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Abstract 

A description of different phases of two dimensional magnetic insu- 
lators is given. 

The first chapters are devoted to the understanding of the symme- 
try breaking mechanism in the semi-classical Neel phases. Order by 
disorder selection is illustrated. All these phases break SU{2) symme- 
try and are gapless phases with ASz = 1 magnon excitations. 

Different gapful quantum phases exist in two dimensions: the Va- 
lence Bond Crystal phases (VBC) which have long range order in 
local S=0 objects (either dimers in the usual Valence Bond accep- 
tion or quadrumers..), but also Resonating Valence Bond Spin Liquids 
(RVBSL), which have no long range order in any local order parameter 
and an absence of susceptibility to any local probe. VBC have gapful 
AS ~ 0, orl excitations, RVBSL on the contrary have deconfined spin- 
1/2 excitations. Examples of these two kinds of quantum phases are 
given in chapters 4 and 5. A special class of magnets (on the kagome 
or pyrochlore lattices) has an infinite local degeneracy in the classical 
limit: they give birth in the quantum limit to different behaviors which 
are illustrated and questionned in the last lecture. 
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Chapter 1 

Introduction 



In this first chapter, we rapidly describe the basic knowledge on Heisen- 
berg magnets to set the frame and the notations of the next developments. 
Different excellent text books can be used for a wider and slower introduc- 
tion [niiEi. 

1.1 History 

The first microscopic model for magnetism goes back to Heisenberg when 
he realized in 1928, that the exchange energy between electrons (introduced 
by Dirac and himself to explain the singlet triplet separation of the gaseous 
helium spectrum) was also responsible for ferromagnetism. Heisenberg and 
Dirac have first suggested that exchange of electrons could be written in an 
effective way in spin space through the use of a spin Hamiltonian, which 
reads: 

h{iJ) = Si.Sj (1.1) 

where Sj, Sj are the spin- 1/2 operators of electrons i and j ^. 

This Hamiltonian allows the basic description of the magnetism of insu- 
lators on a lattice (Heisenberg and Van Vleck). In its simplest form, it is 
written: 

n = jY^ Si.Sj (1.3) 

<i,j> 

where the sum < i,j > runs on pairs of next neighbor sites and J measures 
the strength of the effective coupling ( related to the tunnel frequency of a 
pair of electrons on two neighboring sites). 

""^It may be remembered that h{i,j) is directly related to the spin-1/2 permutation 
operator P{i,j) by: 

P(i,j) = i + 2S,.S,. (1.2) 

Interesting from the conceptual point of view, this relationship is also extremely useful for 
computational purposes. 
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The Heisenberg Hamiltonian, as the more complex schemes of interac- 
tions that will be studied in these lectures, are all SU (2) invariant (i.e. 
invariant in a global spin rotation). On a given sample of N spins, TC com- 
mutes with the total spin of the sample Stot- 

[n,slt] = [n,s^j = o. (1.4) 

The eigen-states of TC can thus be characterized by their total energy and 
their total spin Stot- Eigen-states of 7i with different Stot are a priori non 
degenerate. 

1.2 J < 0: the ferromagnet 

If J is < 0, the ground-state of can be readily written as: 

i 

where the ket |i, -|- > indicates that spin at site i is in the eigen-state of S?, 
with eigen-value +h/2 (in the following we will use h= 1) and the tensorial 
product involves each lattice spin. 

It is easy to show that this state (jl.Sp is an eigen-state of ()1.3|) with the 
total energy: 

Eferro = (1.6) 

where N is the number of spins of the sample, and z the coordination num- 
ber of the lattice. It is also an eigen-state of the total spin Stot and its 
z-component S'i^jwith eigen-values N/2{N/2 + 1) and M^^^ = N/2. It min- 
imizes the energy of each bond and is thus the ground-state of l^l.'A^ . This 
state can be written without ambiguity as: 

\F,>=\Stot = N/2, Mlt = N/2> (1.7) 

It is degenerate with the 2N other eigen-values of Stoti running from N/2 — 1 
to -N/2. 

You should notice that in the thermodynamic limit such a degeneracy is 
negligible: the associate entropy of the extensive ground-state is 0{Ln{N)), 
which does not contradict Nernst Theorem. 

On a macroscopic sample, \Fz > describes a system with a macroscopic 
magnetization pointing in the z direction. Using the total spin operator 
this state could be rotated in any direction u defined by the Euler angles 
{6,<p) with respect to the reference frame ^. We thus obtain the quantum 
description of a coherent state with a magnetization in the u direction: 

^Remember that the third Euler angle x measures an overall degree of rotational free- 
dom ("gauge freedom"), that can be put to in this context. 
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Let me underline that this coherent state remains in the ferromagnetic 
ground-state multiphcity and that its total spin is well defined and equal 
to N/2. 

The semi-classical character of such states is embodied in the follow- 
ing property: the quantum overlap of two states pointing in different di- 
rections decreases exponentially with Stot^ that is with the system size N 
[H]. For macroscopic samples, the state of a ferromagnet can thus be de- 
scribed with classical words and concepts. This can be said in another way: 
the macroscopic spin, understood as a quantum observable, obeys quantum 
commutation relationships: 

[Stot^ ^^ot\ — ^^tot (1-9) 

as Stot is proportional to A^, the relative value of the "quantum fluctua- 
tions" (measured by the commutator/ A^^ ) becomes negligible in the ther- 
modynamic limit. 

The selection of a special state as (|1.5|) to describe the ferromagnetic 
ground-state is a "minor (or trivial?) symmetry breaking" of the problem. 
By SU{2) rotations this eigen-state generates all the ground-state multiplic- 
ity, and this multiplicity only ^ . 

1.3 J > 0: Neel antiferromagnet and spin gapped 
Phases 

1.3.1 A few historical markers 

If J is > the ground-state of is absolutely not obvious. In 1932, Neel 
suggested that the description of experiments was consistent with a picture 
of the ground-state as a special arrangement of ferromagnetic sublattices 
with a zero total magnetization. 

Let us examine the simplest case of the Heisenberg problem on the square 
lattice. This lattice may be partitioned in two sublattices A and B with a 
double unit cell. Each spin of the A (resp. B) lattice is exclusively coupled 
to the B (resp. A) lattice . Such a problem is said bipartite ^. In this case 
we usually write Neel's wave-function as: 

\Cl. Neel w.f. >= JJ |ii , + > |j , - > (1.10) 

i£A,j£B 

^Some authors deny the use of the word "symmetry breaking" in that case, where 
the ground-state, does not involve any mixture of eigen-states with different symmetries. 
They are certainly right from the theoretical point of view. In view of the experimental 
possibility in a macroscopic ferromagnet to point a given direction, we nevertheless use 
this expression with an appropriate qualifier. 

*This is a class of problems, for which exact results are available: Marshall (Peierls) 
theorem, and Lieb and coll theorems. See for example ref 0- 
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This state has indeed a zero component of the total spin in the z direction, 
but it is non zero in the xy plane. This Ising state has maximal sublattice 
magnetizations: Sa,b = Ma,b = The Ising state with a zero 

component of the total spin in the u direction, defined by the Euler angles 
Q and (\), is indeed: 

\Cl. Neel w.f.; u >= e^^'°*'^e''^'°'^|C/. Neel w.f. > (1.11) 

In this antiferromagnetic case the idea that it is possible to restore the overall 
symmetry of the problem by rotation of (|1.1U|) and averaging has more far 
fetching consequences that is usually thought! 

In his biography Neel told that he had to face strong skepticism and 
objections specially from C.J. Gorter (colloquium in Leyden at the Kam- 
merlingh Onnes Lab; 1932). It seems that L. Landau equally rapidly dis- 
carded this special variational wave- function with the same objections as C. 
J. Gorter. 

I have not had access to authenticated sources, but the objections were 
probably of two kinds: 

• The Neel state strongly breaks the SU{2) symmetry of the Hamilto- 
nian and cannot be a good candidate to describe an eigen-state, 

• The existence of ferromagnetic sublattices is not proved and elemen- 
tary quantum mechanics seems in strong disagreement with this as- 
sumption. 

1.3.2 SU{2) symmetry breaking of the Neel states. 

As I will explicit below and explain in details in the next chapter, the Neel 
states breaks the SU{2) symmetry of the Heisenberg Hamiltonian (and the 
lattice geometrical symmetries) in a radical way (quite different from the 
ferromagnetic case). 

The classical Neel state is not an eigen-state of the total spin, and as 
such it can only be described as a linear combination of many eigenstates of 

El. 

In order to have an elementary view of this question let us rephrase Neel 
wave- function in simple quantum terms: two ferromagnetic sublattices A 
and B, defined by their total spins Sa.b = N/A are to be coupled in such a 
way that: 

{S\ + SI^)\Cl. Neel w.f. >=0 (1.12) 

We know from elementary spin algebra that there are (A^/2 + 1) SU{2) 
invariant ways to do this: the different states resulting of this coupling can 
be labeled in an unique way by their total spin Stot , which can range (for even 
N) from Stot = to Stot = N/2. They can be written in an unambiguous 
way under the form ^, Stot, Ms >■ In any of these subspaces, one can 
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indeed select the Mg = component of the total spin, thus fulfilling Neel 
prescription. 

Elementary spin algebra leads to the following expression for the classical 
Neel state wave-function: 

\CL Neel ^.f. >= ( ^''^ _% ^ ) \S..Ss,S,Ms > 

(1.13) 

where 5 runs through the N/2 + 1 possible values of Stot, and in general 
for each value of S, Ms runs from —S to +S. Here the selection rule on 
the z components of the spins implies that Ms = 0. In this expression the 

coefficients | fl^^ 1 are known as Wigner "3i" symbols. These 

\ Ma Mb Ms J 

coefficients are the coefficients of the unitary transformation which trans- 
forms the uncoupled sublattice spins Sa,b to the SU{2) invariant coupled 
combinations. The Wigner "3j" symbols can be calculated by elementary 
algebra, they are tabulated in books and in computer libraries. 

Comparison of this antiferromagnetic coherent state H1.13() to the fer- 
romagnetic one (|1.5|1 . (|1.7|) shows explicit qualitative differences: the fer- 
romagnetic state is a state with a definite total spin N/2, whereas H1.13() 
involves components with total spin ranging from to N/2. This shows 
that Neel wave-function can at best be described as a linear combination of 
a large number of eigen-states of Ti.. 

For bipartite lattices a theorem originally due to Hulthen (1938) j^, Mar- 
shall (1955) in] and strengthened by Lieb and Mattis (1962) [H] states that the 
absolute ground-state of the antiferromagnetic Heisenberg (|1.3|) (and of more 
general antiferromagnetic models respecting the bipartition of the lattice) is 
unique and has total spin zero. Moreover the ground-state energies in each 
S sector are ordered accordingly to Stot- 

ysi, > Stot EoiSi,) > EoiStot). (1.14) 

From that point, we might infer that the \Sa = N/4,Sb = 0,0 > 
state would be a good starting point to describe the absolute ground-state 
l^'o > of (|1.3j) . and forget all the other components of the classical Neel 
state ()1.13|) . But in such a point of view, we lose the foundations for the 
semi-classical approaches: a state with total spin does not allow to point 
any direction in spin space. According to Wigner Eckart theorem, the three 
components of the sublattice magnetizations (as the components of any vec- 
tor) are simultaneously and exactly zero in such a state: ^ 

< ^'o|Sa,b|^o > oc < ^'olStotl^'o > = 0. (1.15) 



''P. W. Anderson and many authors have written that this exact property seems para- 
doxical and contradictory with observations and with the semi-classical approaches (either 
the simplest spin wave approach, as well as, the more sophisticated field theoretical ap- 
proaches laying upon a description of the ground-state by a coherent state). The second 
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To answer Gorter and Landau objection, and support Neel picture for 
quantum antiferromagnets, it is thus necessary to show that eigen-states 
with different Stot appear in the exact spectrum as the different SU{2) 
invariant components of the supposed-to-be quantum Neel state and are 
degenerate in the thermodynamic Umit. In such a Umit, a quantum su- 
perposition of these eigen-states embodies the "strong" symmetry breaking 
associated to Neel's scenario^. Such a mechanism has been described in fuh 
length by P. W. Anderson in two books [3 |H1 IHl • 

In this chapter and for the sake of simplicity only bipartite lattices and 
collinear Neel states are studied.'' 

1.3.3 Ferromagnetic sublattices, "quantum fluctuations" and 
dimer pairing 

The second difficulty with Neel's scenario, the existence of ferromagnetic 
sublattices cannot be supported by quantum mechanics without specific 
calculations. In fact the sublattice magnetizations are not good quantum 
numbers, they are decreased and eventually wiped out by "quantum fluctu- 
ations". This point is common knowledge today. An essential stone mark 
to this understanding is the first spin-wave calculations done in 1952 by 
P.W.A. Anderson HI and R. Kubo 113: ^. 

In this approach one clearly sees that the transverse term of the Heisen- 

assumption is theoretically correct, PW. Anderson knew indeed the answer to the paradox 
and I will describe in the next chapter a simple way to reconcile both approaches. The 
second assumption about experimental observations seems more questionable! This too 
will be briefly discussed in chapter 2. 

^This corresponds to the strict definition of a symmetry breaking situation where the 
macroscopic order parameter does not commute with the Hamiltonian. Technically this 
can happen only by a mixing of different Irreducible Representations (IR) of the bro- 
ken symmetry group. (Elementary example of the broken left-right symmetry in a one 
dimensional problem with an Hamiltonian invariant under reflection). 

^Qualitatively the 3-sublattice Neel state on the triangular lattice has the same prop- 
erties as the collinear state [lOlllllIT^ with the minor difference that the SU(2) invariant 
components of the 3-sublattice Neel states originate from the coupling of three macro- 
scopic spins of length A'^/6. The ground-state multiplicity is thus somewhat larger, and of 
dimension 0{N'^). In this case the demonstration of Lieb-Mattis theorem on the quantum 
ordering of the ground-states energy in each S sector fails: the positive sign property of the 
ground-state wave-function (Marshall property) is no more true. Nevertheless, empirically 
we have observed that the ordering property H1.14^ was realized in exact spectra of most 
systems for large enough sizes: the only restrictions come from systems with competing 
interactions, very near a quantum critical transition to a ferromagnetic state, where we 
have sometimes observed some violations of relation H1.14|l for large S. 

*Even if you are quite familiar with the modern formalism of spin-waves, this paper 
develops a global physical understanding of the subject, and remains an impressive piece 
of work. The conclusion of the 1952 paper of Anderson also describes (in an elusive way) 
the hint toward the solution of the symmetry breaking problem. 
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berg Hamiltonian: 



Si.Sj\i, + > |j, - > = 5f \i, + >\j,-> + 



+ 1/2 



S^ST + SrSl 



\h + > \ > 



= -\ |i,+ >|j,->+^ \i,->\j,+ > (1.16) 

induces spin-flips, decreasing the sublattice magnetizations. These low en- 
ergy excitations (spin-waves) can be described as quantum oscillators: they 
have zero point quantum fluctuations, which renormalize and stabilize the 
Ising energy and decrease the sublattice magnetization. This spin-wave cal- 
culation lays on an ^ expansion and its validity for spins— 1/2 has often be 
questioned. It appears to be qualitatively valid when compared to exact 
results (when they exist) or to more sophisticated numerical work (see Ta- 
ble 1), I will explain in the next chapter the physical reason of this "good" 
behavior. 

To my knowledge exact results exist for 1-dimensional systems (Bethe 
problem, Majumdar -Gosh Ji — J2 problem) where they predict the absence 
of Neel long range order (and algebraic spin-spin decaying correlations in 
the first case, exponentially decaying ones in the second case). For larger 
lattice dimensionality, only the case of the cubic lattice has been shown 
to be Neel ordered |[I4j . On the other hand, the Mermin- Wagner theorem 
precludes existence of Neel long range order (NLRO) at T 7^ for lattices 
with dimension d < 3. This theorem does not give any indications for the 
T = behavior of 2-dimensional magnets which are the central point of these 
lectures. (Rigorous proofs of order exist for spin 1 and larger [T^ [THIIT7] .) 

Neel order versus dimer pairing: naive approach and numerical 
results. 

The classical Neel wave- function (|1.10|) , is a variational solution with an 
energy per bond (— ^). Whereas the quantum ground-state of Hl.l|l is: 

>= > 1^'- > -K'- > U'+ >] (1-17) 

with the energy — | J. This state that we will call either a singlet state or a 
dimer realizes a very important stabilization of a pair of spins (if compared 
to the classical state) but it does not allow to point any direction in spin 
space (it is a state with a total spin zero). At this microscopic scale quantum 
mechanics in its radicalism does not favor the idea of an SU{2) symmetry 
breaking. The controversy on the existence of Neel long range order, specif- 
ically in frustrated (= triangular) geometry or with competing interactions 
has been a long lasting debate opened by P.W.A. Anderson and P. Fazekas 
|18| I19j and fueled again in 1987 with the discovery of High Temperature 
Superconductors in cuprates. 
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Table 1.1: Quantum energy per bond and sublattice magnetization in the 
ground-state of the spin-1/2 Heisenberg Hamiltonian on various simple cells 
and lattices. The sq-hex-dod. is a bipartite lattice formed with squares, 
hexagons and dodecagons. 

When looking in a simple-minded way at a lattice of coordination number 
z, the energy balance between the classical Ising-like Neel state and the 
quantum dimer covering is not so clear. The classical Neel state has an 
energy 

Ed = Y^'-^cosie) (1.18) 

(where 6 is the angle between sublattice magnetizations) to be compared to 
the quantum energy of a dimer covering 

Edim = -y3^. (1.19) 

This simple approach predicts Neel order on the square lattice, it is incon- 
clusive for the hexagonal lattice, or the triangular lattice (which have Neel 
long range order) and it predicts that the Heisenberg model on the kagome 
lattice is disordered (which is correct, see Table [TTT|) . 
Indeed this approach is naive in both limits. 

In the classical limit we have neglected the "fluctuation effects" gen- 
erated by the transverse coupling: these fluctuations effectively contribute 
noticeably to the stabilization of the ground-state of ordered systems (see 
Table irri). 

In the quantum disordered limit, the dimer covering solutions do not 
take into account the resonances between different non orthogonal coverings 
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Table 1.2: The four 2-dimensional phases described in these lectures. 



which are very numerous and are an essential concept for understanding the 
Resonating Valence Bond Spin Liquids (concept introduced in the present 
context by P. W. Anderson in 1973 and named in honor to Linus Pauling). 
The existence of this second kind of phases remained speculative until the 
end of the nineties. We now think that these different scenarios can be real- 
ized in two dimensional spin- 1/2 quantum antiferromagnets (see Table \T7^ . 

In the first part of these lectures, I will try to extract the generic features 
of the Quantum Neel phase in a fully quantum SU{2) invariant framework. 
In so doing I hope to be able to convince you that the symmetry break- 
ing mechanism implemented in the Neel state could be understood from a 
completely quantum and rather simple approach. 

In the second part of the lectures we will discuss the new quantum phases 
where the ground-state does not break SU{2) symmetry and has no long 
range order in spin-spin correlations. We will see that at least two or three 
different phases with these general properties have been exhibited in real- 
istic spin models. The differences between these quantum phases depend 
on the pattern of dimer-dimer correlations: either they display long range 
order and the system is a Valence Bond Crystal or any correlation functions 
are short ranged and it is a liquid (Resonating Valence Bond Liquid). We 
will try to describe the generic properties of their excitations, discuss some 
experimental prescriptions and recent results. 

1.4 Miscellaneous remarks on the use of the words 
"quantum fluctuations", "quantum disorder" 

In antiferromagnets the word "quantum fiuctuations" is often used with 
different acceptions, depending on the context. 

When we say in the spin-wave approach of the antiferromagnet that 
the sublattice magnetization can be wiped out by "quantum fluctuations" 
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let be conscious that it is a model dependent concept! In that case these 
"quantum fluctuations" do not describe a real microscopic, time-dependent 
mechanism: it is just a way to describe a renormalization (we might say a 
dressing) of the Ising-like states. 

On the other hand, when we say in the RVB spin liquid state, that the 
system can fluctuate between different dimer covering conflgurations this 
correspond to true excitations of the system which may be gapped or not. 

Third, when neutronists say that they measured longitudinal or trans- 
verse spin fluctuations, they use the word in its strictest acceptance! The 
root mean square fluctuations of the sublattice magnetization is deflned as 

< ^o\SW^o > . It has the value ^J + 'i^), in the classical Ising-like 
Neel state TTU^ . and is stiU of order 0{N) in any quantum ground-state 
with long range Neel order (as for example, on the square, hexagonal or 
triangular lattice). The same is true of the total spin fluctuations as we 
will understand in the next chapter! And this is observable! The staggered 
susceptibility is the experimental quantity that can be measured experimen- 
tally: it is related to the Fourier transform of the above correlation function. 
It is non zero in NLRO systems and zero in spin gapped ones. 

In fact the fluctuations of the total spin are zero in a "quantum disor- 
dered" system with a spin gap^. But local fluctuations of a configuration of 
spins and dimers in a "quantum disordered" spin liquid can also be observed 
with local probes: as for example muons |25j.... 

Only a few examples of situations that can be uncovered by the loose 
expression "quantum fluctuations"! 



quantum disordered is another awkward expression. The Valence Bond Crystals and 
standard (type I) resonating Valence Bond Spin Liquids are not "disordered systems". 
The degeneracy of their ground-state is lower than the degeneracy of the Neel state and 
they have well defined excitations. Their order is not of the Neel type but they have 
specific order, as we will see in the following lectures. 



13 



Chapter 2 

The semi-classical Neel 
phase: quantum mechanics 
and symmetry breakings 



In this chapter we want to uncover in a very simple quantum mechanical 
point of view, the nature of the semi-classical phase and the ingredients of the 
SU (2) symmetry breaking. This is grounded in the existence of a "tower" 
of SU{2) invariant states which collapse in the thermodynamic limit in a 
ground-state multiplicity that can be described either in the SU (2) invariant 
basis, or in an (overcomplete) basis of semi-classical coherent Neel states. 
We will do it first in a pedestrian calculus approach and then phrase is again 
in a more basic and conceptual point of view parallelizing the translational 
symmetry breaking of solids. The space symmetry breaking will be quickly 
discussed in this chapter. A more detailed study will be done in the next 
chapter where we analyze the mechanism of "order by disorder" in these 
semi-classical antiferromagnets. 

2.1 Calculus approach 

Let us consider the Heisenberg problem on a lattice of N sites with periodic 
boundary conditions. It is interesting to look first to an exactly solvable 
model, that emerges easily from the expression of the Heisenberg Hamilto- 
nian ( Eq. 11.3(1 in reciprocal space: 




(2.1) 



fce-BZ 



In this expression: 




(2.2) 
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where Rj is the coordinate of spin i, N the (even) number of lattice sites 
and k runs on the reciprocal points of the lattice in the first Brillouin zone 
(BZ). 7fc is the structure factor of the lattice: 

7fc = ^ X] cos(k.ei) (2.3) 

i=l,2 

with ej, (j = 1,2), the unit vectors generating the lattice. On this lattice 
the Neel state is invariant by 2-step translations associated to wave-vectors 
k = (0,0) and k = ko = (vr,7r). Let us select these special components in 
the Heisenberg Hamiltonian and rewrite it as: 

n = no + V (2.4) 

with 

Wo = 2J(S2-Sk,.S„k„) (2.5) 

V = 2J 7fcSfc.S_fc (2.6) 

where BZ* is to be understood as the first Brillouin zone minus the k = 
and ko points. 

Simple algebra leads to: 



Ho — 2J(Sq - Skf,.S_ko) 

= ^(SL-Si-S|), (2.7) 

where Stot is the total spin of the sample and Sa,b the total spins of the 
A, B sublattices. 

You might recognize in Tio the toy model used by Lieb and Mattis in 
the demonstration of the ordering theorem jH): it describes a problem with 
constant long range interactions between spins on different sublattices and 
no interactions between spins on the same sublattice. This model can be 
solved exactly. ^ 

^The same kind of toy model can be introduced in the problem of the 3-sublattice Neel 
state on a triangular lattice: in that last case it involves the Fourier components of the 
spins at the three soft points (which are the center and the two non equivalent corners of 
the Brillouin zone) and reads: 

= l^(SL-Si-Sl-S^c) (2.8) 

where Sa,s,c are the total spins of the A, B, C sublattices. Such a model allows the 
same developments as those done below except indeed the comments on the Lieb-Mattis 
ordering theorem 
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2.1.1 The Ising-like Neel state in an SU{2) invariant model 

Hamiltonian Hq (Eq. E^) is an SU{2) invariant Hamiltonian which com- 
mutes with S^oj and Sf^^. It also commutes with S^^, which in this model 
are conservative quantities (good quantum numbers). All these observables 
commute two by two and with TLq . Eigen-states of S^^j, Sf^j, S^, are 
also eigen-states of ()2.7() . with eigen- values: 

4 T 

EiS, Sa,Sb) = ^ [SiS + 1) - SAiSA + 1) - Sb{Sb + 1)] (2.9) 

The quantum numbers for a sample with an (even) number N of sites are: 
. Sa,Sb G [0,1,..,A^/4], 

• For a given set of values of the sublattice magnetizations {Sa, Sb) the 
value of the total spin S £ [\Sa — Sb\, ■■■,Sa + Sb], 

• For a given value S of the total spin of the sample, its z component 
Ms G [-S,-S + 1,....,S -1,S]. 

The ground-state in each S sector Eq{S) is obtained for the maxi- 
mum sublattice magnetization N/i. The energies of these low energy states 
obey the following relation: 

T 4T 
Eo{S) = --(AT + 4) + — [S{S + 1)] (2.10) 

2 

The eigen-states associated to these eigen-values are the SU (2) invariant 
components of the Ising-like Neel state \Sa = N/A,Sb = N/4:,S ,Ms > 
introduced in section H ..S.2I (Eq . II . 13)1 . These eigen-states have four essential 
properties: 

• their number and their spatial symmetries are uniquely defined by the 
coupling of the sublattice magnetizations (these are exact necessary 
requirements) , 

• their sublattice magnetization is N/4, 

• they collapse to the absolute ground-state as 0{j^). 

These levels form a set that has been described by Anderson as a "tower" 
of states [Zl El Ej: we have called them in our original paper QDJS (for 
quasi degenerate joint states). In this lectures we will refer to this set as the 
Anderson tower. 

^This is a special illustration of Eq. I|1.14|l . Let us remark that this ordering property 
is shared by the toy model H2.8|l associated to the 3-sublattice Neel order on the triangular 
lattice. 
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Figure 2.1: Typical spectrum of a finite size collinear Ising magnet. Tlie 
tower of eigen- levels joined by the continuous line and noted |0) is the An- 
derson tower of states needed to form a symmetry breaking Ising ordered 
ground-state (Eg. 12. Ill)) : such a state is non stationary on a finite size sample. 
The second set |1) (dashed line) is associated with the lowest excitations, 
which are highly degenerate and non dispersive. 



On a finite size lattice the classical Neel state ()1.13|) is a non stationary 
state of Ho (eq. 12.7(1 . But, its precession rate decreases as with the 

system size and becomes infinitely slow in the thermodynamic limit. 

The coherent Neel states described by Ea. ()l.ll|) . form an (overcomplete) 
basis of this ground-state multiplicity. The present study of their SU{2) in- 
variant representation shows that the multiplicity of this subspace is 0{N") 
where a is the number of sublattices of the classical Neel state|lHll2j. This 
gives a non extensive entropy of the ground-state at T = in agreement 
with Nernst theorem. 

Excitations 

In this model an excited state is obtained by flipping a single spin of a 
sublattice. From equation ()2.9() one sees that these excitations are localized 
and have an energy: 



rpexc r) T 

^Ising ~ '^•J 



^^4(S + 1) 



N 



f2.111 



For any size these excitations are gapful and 0(J). 
Conclusion 

■^0 describes an Ising magnet in an SU{2) invariant framework: its spec- 
trum has the very simple structure schematized in Fig. 12.11 In the thermo- 
dynamic limit this magnet can be described either in an SU{2) invariant 
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language with the help of the \Sa, Sb, S, Ms > states or with the coherent 
semi-classical Neel states described in Ea. (|l.ll() . The two basis are con- 
nected by exact transformation laws, Eq. (|1.13|) and its inverse: 

\N/4,N/4:,S,Ms >= (23 + 1) j dTDl{^,e)\Cl. Neel w.f .■,u >, (2.12) 

where the differential integration volume reads dr = ■^d(j)d{cosd), where 
(j) £ [0, 27r], 9 £ [0,7r] and Ds is the rotation matrix in the S subspace. In 
the thermodynamic limit, the symmetry breaking point of view is as valid 
as the SU (2) invariant approach. 



2.1.2 "Quantum fluctuations" in the Heisenberg model 

Modification of this picture in an Heisenberg magnet with next neighbor 
exchange comes from the effect of the perturbation term V described in 
Eq. (|2.6|) . V does not commute with and S^: at first order in perturba- 
tion each component of V couples the ground-state of TCq in each S sector 
(which are characterized by maximum uniform sublattice magnetizations) to 
states where the sublattice magnetization is decreased by one unit in some 
modulated way. The analytical treatment of this perturbation is uneasy in 
the SU{2) invariant formalism, but indeed we recognize all the concepts at 
the basis of the usual algebraic spin wave approach: that is renormalization 
of the ground-state energy and of the sublattice magnetization by the zero 
point quantum fluctuations of the spin waves excitations. 
If the structure of the tower of states i.e.: 

• number and spatial symmetries of the states in each S sector, 

• existence in each of these states of a macroscopic sublattice magneti- 
zation (i.e. l^^l > oc 0{N)), 

• scaling as 0{ ■^•^(^+^) ) with respect to the ground-state, 

resists to this renormalization, the nature of the ground-state multiplicity in 
the thermodynamic limit gives a new foundation to the spin-wave symmetry 
breaking point of view The quantum Neel wave function thus emerges 
from the classical picture fEa. ll.ll?)) by the renormalization of the eigenstates 
\Sa, Sb, S, Ms > of Ti.0 under the action of V. We can then write the 
quantum Neel wave-function as: 

IQu. NM ../. >= E ^ ( I X i ) IS:^>0 (2.13) 

b,Ms 

^As an example, all these criteria have been thoroughly checked in the Neel ordered 
phase of the Heisenberg model on the triangular lattice [Tinill| . 
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where the kets |5, Ms >q are now the exact low lying states of the Anderson 
tower oin fEa. O|l . 

2.1.3 The spin- wave algebraic approach 

In order to gather all the material needed for a full understanding of the 
symmetry breaking mechanism in Neel antiferromagnets, let us recall the 
main results of a spin-wave calculation. (For the derivation of the spin-wave 
approach in antiferromagnets, see the above mentioned text-books 0121 E]-) 
Departing from the Ising configuration (Eq ll.lfl|) . the transverse terms 
of the Heisenberg Hamiltonian create = 1 spin flips, which are mobile 
excitations. 

• In an harmonic approximation these excitations are simply described 
as spin-waves, with frequencies: 



2JJI-72 (2.14) 



where 7q is the structure factor of the lattice defined in Eq. (|2.3|) . The 
spin flips excitations are then dispersive, their frequency goes to zero 
when going to the two soft points k = 0, ko. Around these points the 
dispersion law is linear in k (resp. (k — ko)). 

The zero point energy of these excitations (which are oscillator- like) 
renormalizes the Ising classical energy of the ground-state (|1.18|) . To 
first order, this spin wave calculation gives the ground-state energy of 
the Heisenberg Hamiltonian on the square lattice as: 

2 4 ^2 

qG-BZ* 

These " quantum fluctuations" also renormalize the sublattice mag- 
netization. Let us define the order parameter m in the ground-state 
|0 > of this symmetry breaking representation by: 

m=^<Q\S\\Q> . (2.16) 



The first order spin-wave calculation leads to: 

^ (2.17) 



N ^ 



1 



The renormalization of the order parameter is dominated by the fiuc- 
tuations in the low energy modes. The linear asymptotic behavior of 
Ljq around the soft points, implies that the spin-waves correction to 



19 



the order parameter diverges in ID. It gives finite corrections at T = 
on most of the 2-dimensional lattices (square, triangular, hexagonal.. 

• Finite Size Effects: The spin-wave approach allows a direct under- 
standing of the finite size effects in a problem with Neel long range 
order. Let us first remind that on a finite size lattice of linear length 
L, the allowed wave vectors are quantized and of the form This 
introduces a cut-off of the long wave-length fluctuations which is pro- 
gressively relaxed as the size of the sample goes to oo. As is linear 
in q around the soft points, we thus expect that the ground-state en- 
ergy E^~'^ (Eq. I2.15|) and the order parameter m^~^ (Eq. I2.17|) on a 
lattice of finite size L will differ from the L ^ oo limits by factors of 
order 0{j^). This is exactly the result obtained in more sophisticated 
approaches [23 EH EBl OHl EIB EH • 

As we have already underlined the excitations of this model now differ 
from those of Ti.Q: they are itinerant and have acquired dispersion. On 
a finite lattice the energy needed to create the softest excitation is no 
more of order J, but of order ^ oc . 

2.1.4 Self-consistency of the Neel picture for an Heisenberg 
magnet in an SU{2) invariant picture: spectrum and 
finite size effects. 

If the structure of the tower of states is essentially preserved by the quan- 
tum fluctuations due to V, the semi-classical picture of coherent states is 
preserved (see subsection I2.1.2|) , the spin- wave approach is a reasonable one 
and the essential results of this approach should appear in the full spectra 
of Ea. (|1.3l) . Beyond the criteria already described to support the SU{2) 
symmetry breaking, the following size effects should be present: 

• The energy per site of the states of the low lying Anderson tower 
should converge to the thermodynamic limit with a leading correction 
term going as oc 

• The sublattice magnetization ^J~<~\S^\> in each of these states should 
remain 0{N), with a leading term to the finite size corrections of order 

• The low lying softest excitations with wave-vector ^ should be de- 
scribed by a second tower of states issued from the tower of excited 
states of the Ising model with one spin- flip (Eq. I2.11j) . But contrary 

■^The exceptions: the checker-board and the kagome lattice will be studied in a forth- 
coming chapter. 
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Figure 2.2: Typical spectrum of a finite size collinear antiferromagnet with 
Neel order. The tower of eigen-Ievels joined by the continuous line and noted 
|0) is the Anderson tower of states needed to form a symmetry breaking Neel 
ordered ground-state (Eq. 12.13(1 : such a state is non stationary on a finite 
size sample. The second set |1) (dashed line) is associated with the lowest 
magnon. 



to the Ising model, these states are now dispersive and the lowest 
excitation is now distant from the ground-state tower of states by an 
energy of the order of ^: it is the Goldstone mode of the broken SU (2) 
symmetry. 

Some of these properties are summarized in the supposed-to-be spectrum 
of a Neel antiferromagnet described in Fig. 12.21 . This is to be compared 
to an exact spectrum of the Heisenberg Hamiltonian on a square lattice 
(Fig. ESI) '3^ or on an hexagonal lattice fFigs. ITU ^UX ) W' . 

This global understanding of the spectra of finite size samples of anti- 
ferromagnets is a very useful tool to analyze exact spectra of spin models 
that can be obtained with present computer facilities ^. It seems that it may 
equally help to understand the time behavior of nano-scale antiferromagnets 
as ferritin j34j . 

^Historically the first authors to have looked for the Anderson tower of states were 
probably A. Siito and P. Fazekas in 1977 |33| . and with the modern computational facilities 
M. Gross, E. Sanchez- Velasco and E. Siggia |^l?7) . 
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Figure 2.3: Antiferromagnetic Heisenberg model on the square lattice: 
eigen-energies vs eigen- values of S^. The dashed- line is a guide to the eyes 
for the QDJS of the symmetry breaking quantum Neel state (Eq. 12.13(1 . 
The dotted line joins the states associated to the first magnon. There is 
one QDJS for each S (as expected for a coUinear antiferromagnet): they are 
k = states, and k = (vr, tt) states, invariant in C4 rotations. 
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Figure 2.4: Antiferromagnetic Heisenberg model on the honeycomb lattice: 
eigen-energies vs eigen- values of S^. The dashed- line is a guide to the eyes 
for the QDJS. The dotted line joins the states associated to the first magnon. 
There is one QDJS for each 5 (as expected for a collinear antiferromagnet): 
they are k = states, invariant under a 27r/3 rotation around an hexagon 
center, even (odd) under inversion, odd (even) under a reflection with respect 
to an axis going through nearest neighbor hexagon centers for S even (odd) 
(taken from ref. |2()j). 
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Figure 2.5: AF Heisenberg model on the honeycomb lattice, scaling of the 
QDJS with S and iV for = 18, 24, 26, 28, 32 (taken from ref. HU]). 
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Figure 2.6: AF Heisenberg model on the honeycomb lattice, (a) energy 

3 

per site eo versus N~2 (b) spin-gap: The dashed line is a linear fit in 1/A: 
for the sizes of interest the restriction to the leading term of the finite size 
expansion is insufficient. The full line is a fit to eq. [291 I31j : A (A) = 
4^(1 — + where x is the spin susceptibility, c is the spin- 

wave velocity, p the spin stiffness a2d /3 is a number of order one (taken 
from ref. 



2.2 A simple conceptual approach of the transla- 
tional symmetry breaking of a solid 



For a while we exclude any calculations and just rely on very simple and basic 
concepts of condensed matter physics and quantum mechanics to derive the 
"necessary" structure of the spectra of ordered condensed matter in finite 
size samples. For the sake of simplicity, we begin with the problem of the 
solid phase. We successively expose the fundamental classical hypothesis 
underlying the theory of solids. Quantization of this picture enlightens the 
translational symmetry breaking mechanism and finite size effects give a 
new light on the absence of solid order in 1-dimensional physics. 



2.2.1 An essential classical hypothesis 

Let us consider a finite sample of solid with N atoms of individual mass m. 
The Hamiltonian of this piece of solid contains a kinetic energy term and an 
interaction term C/(rj — r^), which essentially depends on distances between 
the atoms, and is translation invariant. Nevertheless any piece of solid in 
nature breaks translational symmetry! 

The first step in the description in classical phase space of the dynamics 
of this object with 2dN degrees of freedom, consists in sorting these variables 
in two sets: 

• the center of mass variables! R-c.o.m 

and Pc.o.m, the dynamics of which 

is a pure kinetic term /C: 

p2 

IC = — — (2.18) 

2Nm ^ ' 



and the 2d{N — 1) internal variables, which obey a dynamic with 
interactions: 

"Hint = 

i€[l,..,N] 



2m 



(2.19) 



Then invoking the inertia principle, the analysis of the problem focuses 
on the Galilean frame, where the center of mass is at rest. In this frame, 
the internal excitations are analyzed in first approximation as modes of 
vibrations: the phonons, which present a dispersion law linear in k for small 
wave vectors k. 

In so doing, an essential dichotomy is introduced between the global 
variable and its dynamics on one hand and the internal excitations on the 
other: this dichotomy is at the basis of the concept of an ordered phase [Hj- 
A technical asymmetry is also introduced in the treatment of the dynamics 
of these two sets of variables: the center of mass dynamics is described in 
a classical framework which explicitly breaks the translation invariance of 
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the total Hamiltonian of the solid /C + 'Hint- On the other hand the internal 
excitations are looked at in a translationally invariant (eventually quantum) 
point of view. This point of view may seem inconsistent in particular when 
looking at a finite sized, eventually small, piece of solid. 

Taking as a definition of the solid phase the essential distinction between 
the global variable and the internal ones, we will show that the technical 
asymmetry in the treatment of these variables can be easily overcome, thus 
explaining both the localization of a piece of solid in real space, and the 
influence of space dimensionality on the definition of this solid. 

2.2.2 Quantization of the classical approach, finite size spec- 
tra, thermodynamic limit and translational symmetry 
breaking 

In order not to break artificially the translational symmetry of the problem 
we consider a solid with periodic boundary conditions. 

If we take for granted that it is legitimate to disconnect the center of 
mass dynamics from the internal excitations we may consider a solid at 
T = with no internal excitations: the vacuum of phonons that we will 
write |0 >. 

The translationally invariant eigen-states of fC are the plane waves with 
wave- vectors k where kx^y^z = nx,y,z^, L is the linear length of the sample 
and nx,y,z non zero integers. Their eigen-values are of the general form: 

iTr- (2.20 

The total energy of the solid in these states is thus of the form: 

^°(') = l^+^- 
where Eg is a constant measuring the zero point energy of the internal 
degrees of freedom. These eigen-states are shown in Fig. 12.71 connected by 
the red continuous line noted |0 >. 

In order to localize the center of mass it is necessary to form a wave- 
packet with eigen-states of /C showing a large distribution of wave-vectors k: 
the largest the k-distribution be, the better the localization of the center of 
mass. Such a wave-packet is non stationary for a finite size, but its evolution 
rate goes to zero as 0{\/N). Localization of the center of mass is thus a 
costless operation in the thermodynamic limit. 

Let us look now to the first excitation of the solid with one phonon of 
wave vector k^m. = 2tt/L. This state can typically be written in a symmetry 
breaking picture as: 

|1 >= exp I ^ik^in.Tj I |0 > (2.22) 
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Figure 2.7: Typical spectrum of a finite size solid. The tower of eigen-levels 
joined by the continuous line and noted |0) is the Anderson tower of states 
needed to form a symmetry breaking vacuum of phonons of the solid: such 
a state is non stationary on a finite size sample. The second set 1 1) (dashed 
line) is associated with the lowest phonon. 

It thus involves a linear superposition of eigenstates of /C + Hint with a dis- 
tribution of wave vectors displaced by k^m. with respect to the distribution 
of the localized ground-state |0 >. This second set of excitations is displayed 
in Fig. 12.71 with a dashed line noted |1 > joining the different eigen-states. 
The softest phonon has an energy proportional to kmin oc oc N~^/'^ 
which should be added to the ground-state energy 1)2. 21(1 giving eigen-states 
with eigen-energies: 

El (k) = —— + Eg + hck, (2.23) 

where c is the sound velocity. Due to the structure of equation ()2.23() the 
line joining the different translation invariant states of this soft phonon is 
parallel to the ground-state line |0 >. This explains the supposed-to-be 
structure of the low lying levels of a finite size solid exhibited in Fig. 12.71 

2.2.3 Thermodynamic limit, stability of the solid and self- 
consistency of the approach 

The consistency of the semi-classical picture implies that the localization 
of the center of mass could be done whatever the degree of excitations of 
phonons: looking to the finite size effects this appears to be the case if 
the dimension of space if larger or equal to 2. In these situations, for large 
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enough sizes there appears two different scales of energy: the Anderson tower 
of states of the ground-state collapses as N^'^ to the absolute ground-state 
whereas the softest phonon collapses on the ground-state only as N~^^^. In 
this limit, the dichotomy between the dynamics of the global variable and 
the internal variables is totally justified. On the other hand in 1 dimension it 
is quantum mechanically inconsistent to separate global degrees of freedom 
from internal ones: these two types of variables having dynamics that cannot 
be disentangled. 

2.3 An analogy: SU{2) symmetry breaking in the 
Neel antiferromagnet 

Let us now develop the analogy between the solid states and the antiferro- 
magnetic ones. 

• The global variables of the solid are Hc.o.m and the conjugate variable 
Pc.o.m- In the collinear antiferromagnetic case the global variables 
of position of the magnet are the two Euler angles {6, (p) allowing to 
point the direction of the sublattice magnetization in spin space. Their 
conjugate variable is the total spin operator S. 

• The free motion of the center of mass is governed by the Hamiltonian 
/C (the quadratic form of this kinetic energy being related to the ho- 
mogeneity of space). By analogy we expect the kinetic energy term 
describing the free precession of the sublattice magnetization to be of 
the form: JCgpin = S^Qj/2wan ^. In such a point of view the constant a 
is just a multiplicative term: we know from other sources (fluctuation 
dissipation theorem or macroscopic approach of the magnet) that this 
is up to a constant the homogeneous spin susceptibility. 

• The eigen-states describing the free precession of the order parameter 
in the vacuum of magnons are thus states with total spin S (ranging 
from to N/2), and eigen-energies: 

Eo{S{S + l)) = ^^^^^+E, (2.24) 

They form the set 1 > of Fig. 12.21 By forming a wave-packet out of 
this set one can localize the direction of the sublattice magnetization 
and break SU (2) symmetry. 

®A three sublattice Neel order has a more complicated order parameter: the three Euler 
angles are needed to localize the 3 sublattice magnetizations: and the macroscopic object 
is no more a rigid rotator as in the case of the collinear Neel order but a (symmetric) 
top. There is in that last case an extra internal spin kinetic energy term and as already 
explained in the previous section the Hibert space of the problem is larger. See ref. 
for example or the quantum mechanical theory of symmetric top molecules. 
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• The discussion of the first excitations above the vacuum of magnon 
completely parallelizes that of the phonons excitations (same disper- 
sion law and same finite size scaling law). The eigen-energies of the 
states embedded in the softest magnon (referred as 1 1 > in Fig. 12. 2() 
are thus of the form: 

E^{S{S + l)) = ^^^^^ + Eg + hc,k^in (2.25) 
2xN 

where Cg is the spin wave velocity. 

• The possibility of a spin rotational symmetry breaking at the thermo- 
dynamic limit is embodied in the finite size behavior of the low lying 
levels of the spectra (Fig. 12. 2|) . In dimension d > 2 the eigen-states 
of the sets |0 > (resp. |1 >) collapse on their S = Smin component 
as 0(A^"i), more rapidly than the decrease in energy of the softest 
magnon which is 0{N-'^/^). In dimension 2 and higher, the SU(2) 
breaking mechanism prevails on the formation of magnon excitations 
justifying the classical approach and the dichotomy between global 
classical variables and internal excitations. 

• These finite size scalings of the Anderson tower of states and of the 
true physical excitations (the magnons) give a new light on the Mermin 
Wagner theorem which denies the existence of Neel long range order 
in 1 dimensional magnets. 

2.4 The coherent quantum mechanical description 
of the Neel state 

At the end of this presentation I hope that Eq. 1)2. 13() now appears as the 
natural quantum mechanical description of a coherent Neel state. And by 
the fact the usual symmetry breaking approach is justified as soon as it gives 
self consistent results (i.e. non zero order parameter). 
The technical answer seems beyond doubt. 

The question is now, do coherent states as those described in Eq. II. IH 
which I rewrite here 

|<3..A..e,».A>=Ei^(«^ _% l)ls3i;>, (2.2.) 
exist in real life? 

I see no mechanism which can lock the difference of phases of the macro- 
scopic number of states \S,Ms >g entering Eq. ()2.26|) to the correct values 
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and it seems that many perturbations could destruct such a coherence, if, 
by an infinitesimal chance, it existed! 

So I will plead that in real life the system may be in any incoherent 

superposition of the degenerate \ S, Ms >q which does not build in spin space 
a given direction to the sublattice magnetization! 

But nobody has to care for it, experiments are not sensitive to the di- 
rection of the sublattice magnetizations but only to correlations functions: 
as the square of the staggered magnetization. This correlation function is 
identical in all the states of the Anderson tower in the Ising model; if it 
survives to quantum fluctuations introduced by V, we expect it to be nearly 

identical in all the \S,Ms >q states at least for total spin up to ~ ^/N 
(above these value of the total spin there might be some difficulties to dis- 
entangle magnons from the Anderson tower of states of the ground-state). 
This has been checked to be true in the Heisenberg model on the triangular 
lattice [Tin . 

As a last remark, the homogeneous spin susceptibility is always domi- 
nated by the largest spin states of the Anderson tower: that is states with 
total spin 0{^/N). Don't forget that a state with total spin \/iV has a 
macroscopic magnetization by site: m = Stot/^ oc that is essentially 
zero in the thermodynamic limit. 

2.5 Space symmetry breaking of the Neel state. 

The Neel state usually breaks some space symmetries of the lattice. 

• In the square lattice case (see Fig. 12. 3p one-step translations are not a 
symmetry operation of the ground-state but the point group is unbro- 
ken. This appears in the Anderson tower of states of Fig. 12.31 where 
the Irreducible Representations (IR) of the QDJS have alternatively 
wave- vector k = (0, 0) or k = (tt, tt) (depending on the parity of the 
total spin), but are trivial IR of the point group. 

• On the hexagonal lattice, which is not a Bravais lattice, the situation 
is somewhat different (see Fig. 12. 5|) : the collinear Neel order does 
not break either the translation group, nor C3, the group of 3- fold 
rotations (noted 7^27r). Only the trivial representation of these two 
groups appears in the QDJS (see Fig. 12. But both the inversion 
group (C2, symmetry operation 7?.^) and the reflection with respect 
to an axis joining the center of the hexagons (a) are broken: these 
symmetry breakings appear in the Anderson tower where there is in 
the QDJS an alternation of even and odd IRs of these two groups. 

Determination of the space symmetries of each S components of the An- 
derson tower can be done exactly using symmetry arguments: the space 



30 



symmetries of each \S, Ms >q depend on S, on the shape and total number 
of spins of the sample 1351 1121 ?] . In the following chapter we will give 
an example of such a determination for the Ji — J2 model on the triangular 
lattice. In a given range of parameters 1/8 < J1/J2 < 1, there is a compe- 
tition between different orders and selection by quantum fluctuations of the 
more symmetric one. The study of this example will show the strength of 
the symmetry analysis and the exact nature of this phenomenon of "order 
by disorder" . 
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Chapter 3 

"Order by disorder 



3.1 Some history 

The concept of "order by disorder" was introduced in 1980 by Villain and co- 
workers|57| in the study of a frustrated Ising model on the square lattice. In 
this model the next neighbor couplings along all the rows are ferromagnetic 
as well as those on the odd columns (named A in the following). The cou- 
plings on the even columns (named B) are antiferromagnetic. It is assumed 
that 

< \Jab\ < Jbb < \Jaa\- (3.1) 

The ground-states of this model have A columns (resp B) ferromagnetically 
(resp. antiferromagnetically) ordered. For a system with a number of sites 

= [mod 4], the degeneracy of this ground-state is 2^^, its entropy per 
spin = -^Ln2 is negligible in the thermodynamic limit. At T = the 
ground-state has no average magnetization and is disordered. The picture 
changes when thermal fluctuations are introduced: it is readily seen that 
a B chain sandwiched between two A chains with parallel spins has lower 
excitations than a B chain sandwiched between two A chains with anti- 
parallel spins. This gives a larger Boltzmann weight to the ferrimagnetically 
ordered system. Villain and co-workers have been able to show exactly that 
the system is indeed ferrimagnetic at low T. They were equally able to 
show that site dilution (introducing non magnetic species) was in a certain 
domain of composition and temperature able to select the same ordered 
pattern, whence the name of "order by disorder" . 

During the nineties several authors have studied a somewhat less drastic 
problem in the classical or quantum Heisenberg model : it is the selection of 
a special kind of long range order among a larger family of ordered solutions 
classically degenerate at T=0 j^, 39, 40. 41 ^.42.. 43, 44_. In the classical mod- 
els, the selection of the simplest ordered structure by thermal fluctuations 
, is due to a larger density of low lying excitations around these solutions, 
whence an increased Boltzmann weight of the corresponding regions and a 
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thermal (entropic) selection of order. 

This same property of the density of low lying excitations can also explain 
a selection of specific spin configurations when going from the T = classical 
approach of the Heisenberg model to the semi-classical one. Suppose that 
many classical spin configurations are degenerate in the classical limit, the 
existence of a larger density of excitations around a specific configuration is 
the signature of a weaker restoring force toward this configuration (larger 
well width in phase space). Insofar as the semi-classical spin- wave approach 
is valid, this implies that the zero point quantum energy X^qg^^. ^ of 
Eq. l2.15l is smaller for this solution, which will thus be energetically selected 
by the "quantum" fiuctuations. Both mechanisms (thermal or quantum) lay 
on the properties of the low lying excitations around the classically T = 
degenerate solutions. 

This selection of order is, in most of the cases, less drastic in the con- 
tinuous spin models, than in the original problem of Villain. In most of 
the cases, the degeneracy of the ground-state is less severe than in the Vil- 
lain case. In the Ising domino problem, the degeneracy of the ground-state 
is 2^ and the thermal selection emphasizes 4 ground-states among these 
2^ . In the Heisenberg problem, as we will see below, in most of the cases 
the less ordered solution has a degeneracy of order 0{N°'), with a the num- 
ber of sublattices, whereas the final order selected by quantum fluctuations 
has only a degeneracy 0{N^) with P < a. From that simple point of view 
one can qualitatively state that the selection of order is less drastic than 
in the Villain problem. A special mention should be done of the Heisen- 
berg model on the kagome, checker-board or pyrochlore lattices. In these 
cases, on which we will return at the end of these lectures, the degeneracy 
of the classical ground-state is exponential in N, there is a residual entropy 
per spin at T = and a selection (if any) of some partial order is a more 
difficult issue. 

3.2 "Order by disorder" in the Ji — J2 model on 
the triangular lattice 

The existence of competing interactions is indeed the main cause of classi- 
cal ground-states degeneracy. As a generic example, one can consider the 
so-called Ji — J2 model on a triangular lattice with two competing antifer- 
romagnetic interactions. This Hamiltonian reads: 

n = 2Ji Si.Sj + 2J2 Yl ^i-^k (3.2) 

<i,j> «i,k» 

where Ji and J2 = aJi are positive and the first and second sums run on the 
first and second neighbors, respectively. The classical study of this model 
has been developed by Jolicoeur et al. j42j . They have shown that for small a 
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Figure 3.1: Top: 4-sublattice classical ground state. Spins in the sublattices 
A and B, as well as spins in C and D, make an angle 29. The plane of the 
spins of A and B makes an angle (j) with the plane of the spins of C and 
D. Bottom: the collinear solutions with the three possible arrangements (in 
this case, classical spins in sublattices A and B are antiparallel) . 



(a < 1/8) the ground state corresponds to a three-sublattice Neel order with 
magnetizations at 120° from each other, whereas for 1/8 < a < 1, there is a 
degeneracy between a two-sublattice Neel and a four-sublattice Neel order 
(see Fig. 13. If) . Chubukov and Jolicoeur [IS? and Korshunov [31| have then 
shown that quantum fluctuations (evaluated in a spin wave approach) could, 
like thermal ones, lift this degeneracy of the classical ground states and lead 
to a selection of the collinear state (see Fig. 13. 1|) [IH] . 

The first study of the exact spectrum of Eq. (|3.2j) done by Jolicoeur et 
al. was not incompatible with this conclusion, but was insufficient to yield 
it immediately. I will show now how the study of the degeneracy of the 
Anderson tower allows a direct derivation of this phenomenon. This part of 
the lecture closely follows the paper by Lecheminant et al. j35j . 

As we have done in section (2.1), let us first study the exactly solv- 
able models which display either four-sublattice order or collinear order. 
These models are obtained by extracting from the Heisenberg Hamiltonian 
expressed in terms of the Fourier components of the spin: 



Sk.S_k 7k + -:7 (cos k. (2ui -|- U2) -|- cos k. (ui -|- 2u2) -|- cos k. (u2 - ui)) 
k ^ 

(3.3) 
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where 7k = l/3^^cosk.u^ (u^ are three vectors at 120 degrees from each 
other and connecting a given site to first neighbors), those which describe 
either the 4-sublattice structure or the colhnear ones. 

3.2.1 Symmetry analysis of the Anderson tower of the 4- 
sublattice Neel order. 

The four k vectors which keep the four-sublattice order invariant are k = 
and the three middles of the Brillouin zone boundaries (cahed in the 
following ki, kn and kc)- It is straightforward to write the contribution of 
these Fourier components to 7i in the form: 

^?^o = ^ (-^1 + J2) (S^ - si - S| - - si) , (3.4) 

where S is the total spin operator and are the total spin operator of each 
sublattice. '^TCq, 8^,8^,8^, 8^ and 8|, form a set of commuting observables. 
The eigenstates of '^TLq have the following energies: 

^E{S, Sa, Sb, Sc, Sd) = ^iJi + J2) [S{S +1)-Sa {Sa + 1) 

-Sb {Sb + 1)-Sc (Sc + 1)-Sd (Sd + 1)] 

(3.5) 

where the quantum numbers Sa, Sb, Sc, Sd run from to and the total 
spin results from a coupling of four spins Sa, Sb, Sc, Sd- 

The low lying levels of Eq. 13.51 are obtained for Sa = Sb = Sc = Sd = 
N/8: 

^Eo{S) = -^^4^ {N + 8) + ^ ( Ji + J2) S{S + 1). (3.6) 

These states, which have maximal sublattice magnetizations S'i = S'g = 
•^c ~ '^D ~ ^("1" + rotationally invariant projections of the bare^ 

Neel states with four sublattices. Their total energy collapses to the absolute 
ground-state as and form the Anderson tower of the 4-sublattice Neel 
order (noted {^E} in the following). 

As we will now show, this multiplicity {'^E} can be entirely and uniquely 
described by its symmetry properties under spin rotations and transforma- 
tions of the space group of the lattice. 

Let us begin by the 8U(2) properties of {^E}. These states result from 
the coupling of four identical spins of length There is Ns different 

ways to couple these 4 spins: the degeneracy of each S subspace is thus 
(25" + l)A'^s') where the {2S + 1) factor comes from the magnetic degeneracy 
of each S eigen-state. Ns is readily evaluated by using the decomposition 
of the product of four spin N/8 representations of SU{2) (P^/^) 

{^E} = ® (3.7) 

^We may say Ising-like Neel state, as these states can be deduced from Ising states of 
the four sublattices pointing in the principal directions of a regular tetrahedron. 
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Table 3.1: Character table of the permutation group S4. First line indicates 
classes of permutations. Second line gives an element of the space symmetry 
class corresponding to the class of permutation. These space symmetries 
are: the one step translation t {A ^ C), 7^27r/3 (resp. '^27r/3) three- 
fold rotation around a site of the D (resp. B) sublattice, and a the axial 
symmetry keeping invariant C and D. N^i is the number of elements of each 
class. 

in spin S irreducible representations {T>^)- One obtains: 

Ns = ^-352 + S(iV + l)+2 + y^ for 5 < ^, (3.8) 

Note that this degeneracy depends both on S and N and not only on the 
total spin S as is the case for a two or three-sublattice problem. 

The determination of the space symmetries of these eigenstates allows 
a complete specification of {^E}. 

• The four-sublattice order is invariant in a two-fold rotation: the eigen- 
states of {^-E} belong to the trivial representation of C2. 

• {^-E} forms a representation of 84^, the permutation group of four ele- 
ments. The eigenstates of {^E^ could thus be labeled by the irreducible 
representations (I.R.) of ^4 (see Table IXT|) . 

• Each element of the space group maps onto a permutation of ^4 : one 
step translations onto products of transpositions as {AC){BD), three- 
fold rotations onto circular permutations of three sublattices {ABC) 
and so on. The complete mapping of the space symmetries of the 
four-sublattice order onto the permutations of 5*4 is given in Table ITT] 
together with the character table of 6*4. 

• Each irreducible representation of 84 can thus be characterized in 
terms of its space symmetry properties. As noted above they are all 
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invariant in TZt^. ri,r2,r3 belong to the trivial IR of the translation 
group, characterized by the wave-vector 0, whereas r4 and have 
a wave-vector kj/jk/ or k^. Ti and r2 belong to the trivial I.R. of 
C3, whereas is the 2 dimensional representation of this same group. 
Finally, Fi and F4 are even under axial symmetry, whereas F2 and F5 
are odd. 

• The number of replicas of Fj that should appear for each S is then 
computed in the S, Ms subspace with the help of the trace of the 
permutations of 5*4: 



n 



LY^TT{Ri\s)ximci{i) 



(3.10) 



where Ri is an element of the class / of 5*4, N(,i{l) is the number of 
elements of the group in this class and Xi(0 the character of the class 
I in the I.R. Fj (see Table m|) . The values of the traces for a given 
total spin S are then found as: 



Tr R. 



Tr R, 



Ms=S, 



(3.11) 



Ms=S+l, 



In each Ms subspace of {'^E} , it is straightforward to find the trace 
of the elements of S4: 



TT^iA,B) (CD) 

Tr ^{A,B) 
TtI{A,B,C,D) 



Ms, 



Ms, 



Ms, 



Ms, 



Ms, 



N/8 

E ^t+v+x+y,Ms 
t,v,x,y=-N/8 

N/8 

t,v=-N/8 
N/8 

E ^3t+v,Ms 
t,v=-N/8 

N/8 

S2t+v+x,Ms 

t,v,x=-N/8 
N/8 

t=~N/8 



(3.12) 



where t, v, x, y, are the z-components of the total spin of each sub- 
lattice (constrained to vary between and —N/8) and 6ij denotes 
the Kronecker symbol. Using equations (|3.1()[ 13.111 13.12j) one readily 
obtains the number of occurrences of each Fj for any S subset of {'^E} 
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N = 16 

Q 


UiZo40D(O 




nriiS) 
nrsiS) 


102021101 
001000000 
202120100 
022322110 
021211000 


N = 28 
S 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 


nr^{S) + nr^{S) 


20515342412 1 1 1 
30425252312 1 
076 11 9 12 9 10 663 3 1 1 



Table 3.2: Number of occurrences nr^iS) of each irreducible representation 
Ti with respect to the total spin S. For N = 28, nr^ and as well as nr^ 
and npj have been added because this sample does not present any axial 
symmetry. 

(Table EHJ. Note that this result depends on S and on the size of the 
sample. 

This symmetry analysis completes the determination of the QDJS of 
{^E}. These properties of the Anderson tower are stable under the action 
of the discarded part of the Heisenberg Hamiltonian. If the ordering of 
levels is not destroyed by quantum fluctuations, the associated quantum 

numbers remain good quantum numbers of the low lying levels \S, Ms >q of 
the Ji — J2 model (|3.2|) . We have thus obtained the complete determination 
(all quantum numbers , and all the degeneracies) of the family of low lying 
levels describing the ground-state multiplicity {'^E} of the four-sublattice 
Neel solutions. 

3.2.2 Symmetry analysis of the QDJS of the Anderson tower 
of states of the 2-sublattice colhnear solutions. 

Let us now consider the collinear solutions (Fig. 13. 1|) . They are particular 
solutions of the 4-sublattice case and we will rapidly get through the same 
scheme of analysis, indicating mainly the new points. The two vectors which 
keep the two sublattices invariant are and the middle of one side of the 
Brillouin zone (the vectors ki, ku and kc correspond respectively to the 
collinear solutions (a), (6) and (c) in Fig. 13.1(1 . Extracting a specific set of 
two wave-vectors from Eq. 13.31 we find the following contribution to the 
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total Hamiltonian: 



no 



The corresponding low energy spectrum for = Sp = N/A is: 



J1 + J2 



{N + 8) + -{Ji + J2)S{S + l) 



(3.13) 



(3.14) 



and is degenerate with the four-sublattice low energy spectrum (see Eq. 13. Hj) . 
But here the two-sublattice have maximal spins Sa = = N/4. These new 
solutions arise from the three symmetric couplings of the 4-sublattice spins: 
Sq = Sa + Sb or Sa = Sa + Sc or = Sa + Sd with the symmetric 
counterparts for Sp. These collinear solutions have thus a Z3 degeneracy 
(see Fig. 13. 1() . The representation space is thus the sum of three products 
jyN/A ^ jyN/A^ -g ^ ^j^g^.^ gjj^pg V^/'^{A,B) (g) P^/^(C,D) and 

T>^/^{A, C)(i^'D^/^{B , D) have in common the same (symmetric) irreducible 
representation with a total spin N/2. On an A^-sample, the representation 
space of the ground state of the collinear solution is: 



5=0 



3V 



5=1 



(3.15) 



The degeneracy is thus 3(25 + 1) for all S values except for S = N/2, where 
it is only (25 + 1). 

The space group analysis is identical to the analysis done for the four- 
sublattice order, but the number of occurrences of each I. R. Tj is now 
different , since the space {'^E} is smaller than {^E}. For each S value there 
are only three replicas of arising from the Z3 symmetry ( Eq. 13.151 and 
Fig. 13. 1|) . This allows the direct computation of the traces of the operations 
of S4 in each S subset of {'^E}. Using the coupling rules of two angular 
momenta (and in particular the fact that the S eigen-state resulting from 
the coupling of two integer spins changes sign as (— 1)*^ with the interchange 
of the two parent spins) one obtains (for S 7^ N/2): 



Tr{Id\s) = 3 

TriiA,B){C,D)\s) = l + 2{-lY 
TriiA,B,C)\s) = 
TT{iA,B)\s) = l 
TT{{A,B,C,D)\s) = {-lf 



(3.16) 



for 



Therefore the collinear solution is simply characterized by Fi and T3 
even S, and F4 for odd S, whatever the sample size. 

The symmetries of all states of the tower are now fully determined both 
for the 4-sublattice order {'^E} and for the collinear order {^E}. If the 
quantum Hamiltonian presents one of these kinds of order, the quantum 
fluctuations generated by the discarded part of Ti. should preserve the dy- 
namics and the structure of these low lying subsets. 
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Figure 3.2: Top: complete spectrum of the = 16 periodic sample with 
respect to S^. Bottom: enlargement of the difference between the exact 
spectrum and the energy of the low lying levels of the model Hamiltonians 
f Eg 13.61 or Ea l3.14() . The ground-state multiplicity {^E} is well separated 
from the magnons. 

3.2.3 Exact spectra of the Ji — J2 model on small samples 
and finite size eff'ects: a direct illustration of the phe- 
nomenon of "order by disorder". 

We have determined the low (and high) energy levels of the Ji — J2 Hamil- 
tonian in each I. R. of SU (2) and of the space group of the triangular lattice 
for small periodic samples with A'^ = 12, 16 and = 28. The spectra are 
displayed in Fig. 13.21 and Fig. 13.31 We directly see in the upper parts of 
these figures the set of QDJS ("Anderson tower of the ground-state") well 
separated from the set of levels corresponding to the one magnon excita- 
tions. We have verified that this set has the symmetry properties of the 
above defined {^E} subset. The action of the quantum fluctuations could 
then be read in the lower parts of the figures. As expected, the quantum 
fluctuations lift the degeneracies which are present in the exactly solvable 
model and stabilize the eigenstates with the lower S values. Nevertheless the 
low lying energies per site still group around a line of slope 0[{Ji + J2)/N'^)- 
The number and space symmetries of these levels for each S and A^ value 
are exactly those required by the above analysis of the four-sublattice Neel 
order. 

Moreover, it is already visible on the A" = 16 sample and quite clear on 
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Figure 3.3: Partial spectrum of the = 28 periodic sample. (Same legend 
as for fig l3.2|) . Bottom: the tower of states of the 4-sublattice order {^-E} 
lays under the dashed line. Above appear the first magnons. Above the 
dotted line are represented the first excited homogeneous states. In the 
magnon multiplicity (k / 0,kH,ki or kc), for 5" < 5, only the lowest 5 
states of each I.R. have been computed. 



41 




Figure 3.4: Enlargement of the = 16 and = 28 QDJS. A global con- 
tribution I3Eq{S) is subtracted from the exact spectrum. This contribution 
describes the overall dynamics of the order parameter in this finite sample, 
(3 measures the renormalization of this dynamics by quantum fluctuations 
(see refs. ^2 EEj)- The bars represent eigenstates which belong both to 
and {^-E}. The triangles indicate states which belong to {^i?} but not 
to {^-E}. With increasing sizes, the subset {^i?} is stabilized and separates 
from the pure 4-sublattice order. For N = 28 the two states of {"^E} with 
even S are quasi degenerate and cannot be distinguished at the scale of the 
figure. 

the = 28 sample that a dichotomy appears in this family (see Fig. 13.4)1 . 
The lowest levels of this tower of states appear to be FijFs or r4 represen- 
tations depending on the parity of the total spin. They precisely build the 
family {'^E} of isotropic projections of the collinear solutions (Eq. I3.16|) . 
This strongly suggests that the 4-sublattice order will disappear in the ther- 
modynamic limit and only the collinear order will subsist, as was predicted 
in the spin- wave approach |421 1431 144j . 

3.3 Concluding remarks 

• The symmetry and dynamical analysis of the low lying levels of a 
Hamiltonian likely to exhibit ordered solutions gives rather straight- 
forward answer to the kind of order to be expected. The method is 
rapid, powerful and unbiased. It does not require any a priori sym- 
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metry breaking choice: if a specific order is selected, one should see 
it directly on the exact spectra. Moreover, as it is essentially exact, 
there are no questions relative to the convergence of the expansion as 
in the spin-wave approach. On the other hand, as the sizes amenable 
to computation are limited, there is, in the exact approach, a cut-off of 
the long wavelength fluctuations. Results so obtained should thus be 
examined in the light of a finite size scaling analysis. This work never- 
theless shows that it is not necessary to invoke quantum fluctuations 
with very long wave-lengths to select the collinear order. 

• The selection of "order by disorder" appears in a particular clear light. 
Increasing the sample size, increases the presence of long wave-length 
fluctuations. We see on this example how these long wave-length fluc- 
tuations realize a differential stabilization of the {'^E} subset, favoring 
collinear order and progressively wiping out 4-sublattice order. This 
fully support the spin-wave calculations. This is also a clear illustra- 
tion of the previous comment on the "non drastic" character of this 
phenomenon in this peculiar case. Without fluctuations the system 
has already some order: the role of the quantum fluctuations is just 
to restore a higher degree of symmetry to the ground-state solution. 

• Going along this route one may always speculate if in the thermody- 
namic limit, quantum fluctuations could not completely restore the 
symmetries of the Hamiltonian. Spin-waves calculations, so long they 
are consistent at small sizes with exact diagonalizations and self-consistent 
when going to the thermodynamic limit are credible. This compari- 
son is always useful and relevant: in the Ji — J2 — J3 model on the 
hexagonal lattice \20, ?], there are regions of parameter space where 
finite-size exact diagonalizations give Neel Long Range Order whereas 
spin-waves in the thermodynamic limit indicate an absence of sub- 
lattice magnetization! We have verified in each of these cases that 

at small sizes the semi-classical solution in the spin-wave approach 
was equally robust and was only destroyed by very long wave-length 
fluctuations. 

• On the other hand, in the situations where we claim an absence of 
Neel Long Range Order and a more exotic phase (see next chapters 
on Valence Bond Crystals and Resonating Valence Bond Spin Liquids) 
the Anderson tower of states is absent even on the smallest sizes. In 
such a case it is quite clear that the spin-wave approach should be 
discarded for spin-1/2. 
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Chapter 4 

Valence Bond Crystals 



4.1 Introduction 

In our quest of exotic quantum ground-states, we will now describe some 
examples where the semi-classical Neel order is not the ground-state of the 
problem and SU{2) symmetry is not broken. 

In this chapter we will concentrate on solutions where there is long 
range order in the dimer coverings: we call these phases Valence 
Bond Crystals (in the following noted VBC). 

Such solutions are well known in 1-dimensional problems as for example 
in the A.F. Ji — J2 model: 

n = JiYl Si.Sj + J2 s-^i (^-1) 

<ij> 

where the first (resp. second) sums run on first (resp. second) neighbors. 
In 1-d, for J2/J1 > 0.24, the ground-state is dimerized and there is a gap to 
the first excitations: this is the simplest case of a VBC. 
What is the situation in 2-d? 

In a classical approach, the ground-state of Eq. 1)4. 1() on a square lattice 
has a soft mode at (vr, vr) for J2/J1 < 0.5. At J2/J1 = 0.5, the ('/r,7r) order 
is degenerate with 4-sublattice order and collinear (vr, 0) or (0, vr) order. For 
J2/J1 > 0.5, quantum fluctuations select the collinear (vr,0) or (0,vr) order 
by the phenomenon of "order by disorder" (see Fig. 14. 1() . 

In our naive approach of chapter 1, comparing the energies of classical 
Neel solutions to dimer covering ones, we would conclude that dimer covering 
solutions and VBC are more stable than any classical Neel order in a large 
range of parameters around J2/J1 = 0.5 (Fig. 14.1)) . 

In fact "quantum fluctuations" stabilize the Neel states and the window 
for an exotic phase is smaller than indicated in Fig. 14.11 The nature of the 
quantum phase on the square lattice at J2/J1 = 0.5 is still debated 1171 
I48| I49 | 1^ I51j. A columnar VBC has been identified in the same model on 
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Jl/Jl 



Dimers 



Selection of order by disorder 
&; Restoration of symmetry by fluct. 




columnar VBC staggered VBC 

Figure 4.1: Schematization of different variational solutions of the Ji — J2 
model described in the introduction of this chapter 
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Figure 4.2: The checkerboard lattice: the spins sit at the vertices shown by 
bullets, all couplings are identical, ui, U2 are the unit vectors of the Bravais 
lattice. 

the honeycomb lattice for Jij J\ ~ 0.4 (see ref. |2Uj and refs. therein). For 
a pedagogical illustration we will move to a more clear-cut example: the 
Heisenberg model on the checkerboard lattice |52| ?] (noted in the following 
HCKB). 

4.2 The Heisenberg model on the checker-board 
lattice: an example of a Valence Bond Crystal 

The checker-board lattice is made of corner sharing tetrahedrons, with all 
bonds equal: this a 2-dimensional slice of a pyrochlore lattice. The under- 
lying Bravais lattice is a square lattice and there are two spins per unit cell 

(Fig. 1321). 

4.2.1 Classical ground-states 

The Heisenberg Hamiltonian on such a lattice is highly degenerate in the 
classical limit. Due to the special form of the lattice this Hamiltonian can 
be rewritten as the sum of the square of the total spin of corner sharing 
units a : 

n = j - ^ E - (4.2) 

bonds a units 

A classical ground-state is obtained whenever Va = 0. Such ground- 
states have a continuous local degeneracy and an energy —{NJ)/4. This is 
much higher than the dimer covering energy, which is — (3iVJ)/8. As we will 
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Figure 4.3: Gap between the absolute ground-state and the first S=l exci- 
tation of the HCKB model versus sample sizes. 

see below, there is no memory of these classical solutions in the quantum 
ground-states and low lying excitations of this model. 

4.2.2 The Quantum HCKB model: Spin Gap 

As we have seen in chapter 2, the first characteristic of the semi-classical 
Neel like solution is the existence of the Anderson tower of states which 
collapse to the ground-state as 0{1/N) and the absence of spin gap in the 
thermodynamic limit (see for example Fig. 12. 5|) . 

The first salient feature of the Heisenberg model on the checker-board 
lattice is the existence of a large spin gap, which shows no tendency of going 
to zero at the thermodynamic limit (compare Fig. 14.31 with Fig. 12. 5|) . This 
indicates that the ground-state does not break the SU{2) symmetry of the 
Hamiltonian, and as a corollary we expect that the spin-spin correlations 
decrease to zero at large distance (which seems well verified, see Table IV 
of ref. [?]). 

4.2.3 Degeneracy of the ground-state and space symmetry 
breaking in the thermodynamic hmit 

The low lying levels of the spectra of the HCKB model in the singlet space 
are displayed in Fig. 14.41 
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Figure 4.4: Gaps to the first (open up triangles) and second (black squares) 
level of the singlet sector. For the studied samples these two "excited" 
singlet levels are in the singlet-triplet gap (See Fig. I4.3|l . 



In this figure, one reads that the first excited singlet state very plausibly 
collapses to the absolute ground-state, whereas a finite gap to the third S=0 
level (perhaps smaller than the spin gap) build on with sample size. This 
pleads in favor of a 2-fold degeneracy of the absolute ground-state in the 
thermodynamic limit. 

The absolute ground-state is in the trivial representation of the lattice 
symmetry group. Its wave function is invariant in any translation and in 
any operation of D4: group of the 7r/2 rotations around point O (or any 
equivalent point of the Bravais lattice) and axial symmetries with respect 
to axes ui and U2 (see Fig. 14. 2|) . The excited state which collapses on it in 
the thermodynamic limit has a wave vector (vr, vr) (its wave function takes 
a (-1) factor in one-step translations along ui or U2), and it is odd under 
7r/2 rotations and axial symmetries. In the thermodynamic limit the 2-fold 
degenerate ground-state can thus exhibit a spontaneous symmetry breaking 
with a doubling of the unit cell. 

Such a restricted symmetry breaking does not allow a columnar or stag- 
gered configuration of dimers: both of these states have at least a 4- fold 
degeneracy fFig. 14.^ . The simplest Valence Bond Crystals that allow the 
above-mentioned symmetry breaking are described by pure product wave- 
functions of 4-spin 8=0 plaquettes. 

This family includes eight different configurations: 

• The singlet plaquettes may sit either on the squares with crossed links 
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Figure 4.5: Columnar and staggered configuration of dimers (fat links) on 
the clieckerboard lattice: such symmetry breaking configurations are 4-fold 
degenerate in the thermodynamic limit. 




5 T ST 



A configuration B configuration 



Figure 4.6: S=0 4-spin plaquette valence-bond crystals on the checkerboard 
lattice: fat links indicate 4 spins involved in a singlet. 

or on the void squares (A and B configurations of Fig. I4.6|l . 

• The translation symmetry breaking configurations may be in two dif- 
ferent locations named >ii(2) (resp -Bi(2))) 

• An S=0 state on a plaquette of four spins sitting on sites (a, (3, 7, 6) 
may be realized either by the symmetric combination of pairs of sin- 
glets: 

>= \a ^ 6 > \ j ^ P > +\a ^ /? > I7 ^ (5 >, (4.3) 
or by the anti-symmetric one: 

IV"" >= |a ^ (5 > I7 ^ /3 > -\a P > \'y ^ 6 > . (4.4) 
where |a — > 7 > is the singlet state on sites a and 7: 

|a^7>= (|aT,7i> -|ai,7T>)/^- (4-5) 
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Tn = B^ + r^B^ 




(-)Pr?r'' 





Table 4.1: Transformation rules of the product wave- functions in the ele- 
mentary operations of the symmetry group (the space group is defined with 
respect to point O and translations ui, U2). The wave- functions of the anti- 
symmetric plaquettes have different symmetries depending on the parity p 
of the number of plaquettes in the sample. 

We can thus define eight different product wave-functions labeled: |^^2) ^ 
and 1-6^2) >■ The transformations of these states under the elementary op- 
erations of the lattice symmetry group are described in the first four lines 
of Table 14. 11 The symmetric (resp. anti-symmetric) linear combinations of 
these states which are irreducible representations of this group are defined in 
the four last lines of the same Table. Comparison of the symmetries of these 
states for different samples with those of the two first levels of the exact 
spectra indicates a symmetry of the HCKB ground-state doublet. 

In the thermodynamic limit the symmetry breaking configuration is thus of 
the B type decorated by the symmetric 4-spin plaquettes > described 
in Eq. Ol 

A simple last remark could be done: the symmetric-plaquette state 
(Eq. 14. 3() can be rewritten as the product of two triplets along the diag- 
onals of the square. This configuration of spins is not energetically optimal 
on the squares with antiferromagnetic crossed links (A configuration) but 
might a priori be favored in B configuration. Reversely the -plaquette can 
be rewritten as the product of two singlets along the diagonals of the square, 
and would eventually be preferred in A configuration. The variational energy 
per spin of the product wave-function of "0 "'^-plaquettes in B configuration 
is E^ar{B~^) = —0.5, whereas the variational energy per spin of the product 
wave-function of V'" -plaquettes in A configuration is E.uar[-A~) = —0.375. 
The exact energy per spin is E^x ~ —0.514 it 0.006. This is a first proof 
that the real system takes advantage of some fluctuations around the pure 
product wave-function to decrease its energy. 
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Figure 4.7: Dimer-dimer correlations in the exact ground-state of the 36 
sample fEq. 14. Hj) . The reference bond is the bond (1,2). Positive (negative) 
correlations are drawn as full (dashed) lines. The thickness of the lines is 
a measure of the strength of the correlation. The diagonal lines show the 
position of the crossed links. 

The study of dimer-dimer correlations (Fig. 14.71 and Table IX2|) : 

C^{l,2;i,j) =4[< S1.S2 S,.Sj >-< S1.S2 >< Si.Sj >] (4.6) 

and 8-spin correlation functions [?] shows long range order in the 4-spin 
plaquettes, but also the dressing of the pure product state Z+ by quantum 
fluctuations (see Table . 

Are those small size computations relevant for the description of the 
thermodynamic limit? The stronger answer is read in Fig. 14.31 & Fig. 14.41 
insofar as the degeneracy of the ground-state and the gaps to the first triplet 
state and the third singlet state remain finite in the thermodynamic limit, 
the Valence Bond Crystal picture (with LRO in plaquettes) will survive to 
quantum fiuctuations. The gaps results (Fig. 14.31 fc Fig. 14. 4j) show that the 
studied samples (except the N = 16) have linear sizes of the order of, or 
larger than the spin-spin correlation length. We thus think that the present 
qualitative conclusions are reliable. 
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Table 4.2: Dimer-dimer correlations C^(l,2;i,j) (Eq. 14. 6|) in the A'^ = 36 
ground-state. The sites l,2,z, j are described in Fig. 14. 7( the i,j points are 
enumerated in the first columns. This correlation has been measured in the 
exact ground-state wave function (second columns) and in the variational Z 
state (third columns). All the values of these correlations between sites of 
Fig. 14.71 can be obtained from this table by a mirror symmetry through the 
bisector of bond (1, 2). 
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4.2.4 Excitations: raw data and qualitative description of 
the first excitations 

Looking to Table Ol and Fig. ESI it appears that the triplet excitations 
are gapped (gap of the order of 0.7) and very weakly dispersive. Singlet 
excitations too are gapped (4th line of Table 14.31 and Fig. 14.9)1 : they are 
much more dispersive than the triplet excitations and less energetic (gap of 
the order of 0.25). 
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Table 4.3: Spectrum of the Heisenberg model on the checker-board lattice. 

n' 

Energy per spin in the ground-state cq and energy gaps - E^f between 
the n'g energy level of the S' spin sector and the ng level of the S sector. 
Second line: spin gap. Third line: gap between the absolute ground-state 
and the first singlet excitation . Fourth line: gap between the second and 
third level in the S = sector. Fifth line: gap between the third level 
in the 5 = sector and the first triplet excitation. Following lines: ni is 
the number of singlet states in the spin gap (including degeneracies). The 
starred column corresponds to a sample which has the extra symmetries of 
the pyrochlore lattice. 

There is a very simple variational description of the triplet excitations: 
let us consider the 4-spin plaquettes B of the ground-state. The S=0 ground- 
state is formed from the coupling of two triplets along the diagonals. There 
are four S=l states on such a plaquette. The lowest S=l excitation simply 
results from the S=l coupling of the two diagonal triplets. The gap to this 
variational excitation is 1. The Bloch waves built on such excitations are 
non dispersive. Up to a renormalization of the gap of the order of 33%, this 
picture appears as a good qualitative description of the true S=l excitations 
of the HCKB model, which are massive, quasi localized excitations with an 
energy gap ~ 0.7. 

The singlet excitations are more intricate. On a B plaquette the first 
S=0 excitation corresponds to the antisymmetric coupling of dimers \ip~ > 
described in Eq. 14.41 Its energy gap to the ground-state is equal to 2. This 
first excitation of the B plaquettes is more energetic than the object built by 
a reorganization of two symmetric ip'^ states on two neighboring B positions. 
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Figure 4.8: Dispersion relations in the triplet sector versus |k|/|ko| with 
ko = (7r,7r). 

The inset shows the correspondence between the colors of the symbols and 
the wave vectors in the Brillouin zone. Only the triplet excitations are drawn 
in this figure. 
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Figure 4.9: Dispersion relation of the singlet excitations of the A'^ = 36 sam- 
ple versus |k|/|ko| with ko = (7r,7r). The horizontal dashed line indicates 
the spin-gap. The inset shows the correspondence between the colors of the 
symbols and the wave vectors in the Brillouin zone. Only the singlet exci- 
tations are drawn in this figure. The first excited level of this figure with 
ko = (tt, tt) (yellow up triangle) is not a true excitation. This level is at the 
thermodynamic limit degenerate with the ground-state and allows the space 
symmetry breaking of the 4-spin plaquette order. 
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More precisely the excitation which promotes the two pairs of spins (a, 6) 
and 7) into triplet states and then couples them in a singlet states has a 
gap 1 with respect to the ground-state. To first order in a strong coupling 
expansion this excitation is non dispersive but it can acquire dispersion 
at higher order. The exact S=0 excitations are thus certainly a bit more 
extended and complex that this first approximation.^ 

Remarks on the number of singlet excitations in the singlet- 
triplet gap 

The spectra of very small samples of the HCKB model j52| lead to con- 
clusions on the number of singlets excitations that should be precised and 
relativized. The above-mentioned authors, and many commentators after- 
ward, argued that this large number of singlets might be reminiscent of the 
exponential degeneracy in the singlet sector of the Heisenberg model on the 
kagome lattice. A precise analysis of the spectra of singlet excitations of 
the HCKB model shows that this analogy is unjustified for the following 
reasons: 

i) The sample set used to extract this conclusion was a mixture of true 
2-d checker-board samples, quasi 1-dimensional tubes with a cross section of 
4 spins and 3-dimensional pyrochlore-like samples [?]. Table ^21 summarizes 
the results for the restricted family of pure 2-dimensional HCKB samples. 
Whereas the low density of singlet states of the Heisenberg model on the 
kagome lattice increases as 1.15^ with the system size, a similar analysis 
for the HCKB model gives the value 0.14''^. Changing the ratio of the 
exponential from a number larger than one to a number smaller than one 
changes indeed completely the picture! In fact the best fit to describe the 
number of singlet in the spin-gap is obtained for a power law fit: with 
7 = 3.96. If these low excitations can be described as modes one would 
effectively expect some power law. 

ii) Last difference between the HCKB model and the Heisenberg model 
on the kagome lattice (HK): in the HK model the continuum of singlets is 
adjacent to the ground-state whereas, in the HCKB model, there is a clear- 
cut gap in the singlet sector between the ground-state multiplicity and the 
first singlet excitations. 

^In view of the strong VBC correlations of the ground-state different authors developed 
strong coupling perturbative studies of the excitations. The work by Brenig et al. |53| 
takes as a departure point the 4-spin S=0 plaquettes on the B positions and treat the 
couplings between the B sites as a perturbation. The limit toward the isotropic point 
of this high order perturbation expansion seems to be rather well behaved but it fails 
to restore the correct symmetry breaking of the ground-state in the S=0 sector. This 
is perhaps not too surprising in view of the above remarks on the first excitations in 
the singlet sector. Berg and collaborators [?] have used a more sophisticated method 
(determination of an effective hamiltonian by a real space renormalisation method called 
CORE), with a different departure point in the S = sector, they found domain walls 
between the two degenerate ground-states as the lowest excitation in the S = sector. 
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This simple analysis of the singlet and triplet excitations of the HCKB 
model supports the idea that the excitations are simple "optical modes" that 
could be observed in Raman, RPE, ESR or NMR spectra. This structure 
is certainly highly reminiscent of those of dimerized or spin-Peierls chains, 
gapped ladders.. A picture consistent with the strong coupling description 
of the ground-state. 

4.2.5 Summary of the generic features of a Valence Bond 
Crystal 

I would like to argue that the main features of Valence Bond Crystals (what- 
ever the dimensionality of space where they are living) are probably: 

• A spin gap, no SU{2) symmetry breaking and short range spin-spin 
correlations, 

• Degeneracy in the thermodynamic limit of the S = ground-state, 
embedding the spontaneous space symmetry breaking of the phase, 

• Long range order in dimer-dimer and/or larger 5 = plaquettes, (that 
can subsist up to finite temperature?) 

• Gapped excitations, in the S = sectors as well as in other S sectors, 
that can be described as modes, more or less dispersive. A strong 
coupling analysis of these modes seems a priori valid but the examples 
worked out up to now, on this model or on the Shastry-Sutherland 
model [SSI innii show that the departure point of the perturbation 
theory should be given special consideration. 

No experimental evidence exists up to now of a pure Valence Bond Crys- 
tal with spontaneous symmetry breaking. But a few 2-dimensional sys- 
tems with a Valence Bond ground-state have been observed experimentally: 

caV409 inziEEiEsinoiiniiiniiMiiMiinsiEiiniiinHiisni and srCu2(B03)2 ^ 

ESlIZIlIZaiZHlElilZnilZHlinElIZZlIZHl- However, in both cases the ground-state 
is non-degenerate because the Hamiltonian has an integer spin in the unit 
cell (4 spins 1/2) and the dimerization does not break any lattice symmetry. 
Prom the theoretical point of view on may argue that the difference between 
the theoretical VBC described in this chapter and these compounds is akin 
to the difference between the dimerized phase of the Ji — J2 model on the 
chain and the spin Peierls compounds. This is a minor difference and from 
the experimental point of view the first excitations of all these models can 
be qualitatively described as "optical modes". The detailed characteristics 
of the two-quasi-particle continua could be a bit more different. 
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4.3 A simple model of VBC with a critical point: 
the hard core quantum dimer model of Rokhsar 
and Kivelson on the square lattice 

Looking for a model with a resonating valence bond ground-state, Rokhsar 
and Kivelson introduced in 1988 a quantum hard core model on the square 
lattice |79j . Their motivation was the description of systems with strongly 
coupled real-space Cooper pairs. At half filling these next-neighbor Cooper 
pairs can be seen as next-neighbor dimers. Pauli principle and Coulomb 
interaction imply that these dimers are hard core dimers. Insofar as the 
spin gap is large enough, it can be speculated that the manifold of low 
energy states is spanned by the linearly independent set of nearest neighbor 
dimer coverings^ . The dynamics of the low lying singlet excitations of this 
model are described by the Hamiltonian TLdimer'- 

H,imer= ^ ["^ (I I I ) I + M + ^ (I I I ) (l I I + I H ) I ) ] 

Plaquette 

(4.7) 

(In their original paper the authors discussed the derivation of this effective 
Hamiltonian from a more realistic Hubbard model.) 

The first term of Eg. 14.71 describes the spatial flip of two parallel dimers 
from horizontal to vertical position and vice-versa, it could also be seen as a 
cyclic permutation of the two dimers around a square: it is a kinetic energy 
term which favors resonances between different configurations of parallel 
dimers (J is always > 0). The second term is a potential energy term likely 
to be repulsive in the original electron model. The ground-state for infinitely 
large is a Valence Bond Crystal, either staggered (for large ^ > 0), or 
columnar (for large -j < 0) . See Fig. 14.31 

Topological structure of the Hilbert space of the QHCD model 
on the square lattice 

The eigenstates of TLdimer can be classified according to their winding 
numbers {Vtx-,^y) across the 2-torus of the square sample with periodic 
boundary conditions. There are many equivalent ways to define these wind- 
ing numbers. Let us follow RK. They draw the transition graph of a dimer 
configuration C relative to a reference configuration Cq (which may be the 
columnar configuration) as the superposition of the dimer coverings of the 
two configurations C and Cq. The dimers in C are directed from one sublat- 
tice to the other and reversely for the dimers of Cq. The transition graph thus 
appears as a graph of oriented loops. The winding number Vlx (resp. ^y) 
measures the net number of loops (clockwise minus counter-clockwise) encir- 
cling the torus in the x (resp y direction) . The Hamiltonian does not couple 

■^It has been shown that such set form a family of non-orthogonal but linearly indepen- 
dent states IHOIIHI]- 
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Rokhsar-Kivelson hard core quantum dimer model 
on the square lattice '88 
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Similar model on the hexagonal lattice; Moessner and Sondhi '01 



subspaces with different winding numbers. These pairs of winding numbers 
define Ns disconnected topological subspaces (where Ng is the number of 
lattice sites). 

For y > J > 0, Hd imer IS positive semi-definite, the ground-state is 
unique and nodeless (Frobenius theorem). Moreover is a lower bond of its 
energy. 

Demons: A lower bond for the ground-state energy is given by a min- 
imization of the Hamiltonian on each plaquette individually. If the given 
plaquette has no parallel dimer (non flippable plaquette), its energy is zero 
and if it has parallel dimers it has a potential energy V and at best a ki- 
netic energy of — J. We can thus write a lower bond energy of the global 
system, which is proportional to the number of flippable plaquettes rifnp, as 
min [0, {V - J)nfiip]. 

Phase Diagram of the RK model 

The staggered configuration is a zero-energy eigenstate of TCdimer ■ 
At the point -j = 1 the model is exactly solvable. 

• The four staggered configurations are zero-energy eigenstates oiTidimer- 
As they saturate the low energy bond, they are the ground-states for 
y > J > 0. They can be classified in two different topological classes in 
which they are the only representatives. They have a zero energy and 
any configuration of other topological subspaces has a larger strictly 
positive energy (at least of order 0{L) in the limit y ^ oo ). 

• At the point y = 1 the model is exactly solvable. There is in each 
topological subspace a ground-state with zero energy. It is the equal 
amplitude superposition of all the configurations of that sector, i) 
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Simple computation shows that these states are zero-energy states of 
'Hdimer ■ H) Since all off-diagonal elements are non-positive the ground- 
state is unique and nodeless (Frobenius property, Marshall theorem). 
The equal amplitude states are thus the unique ground-states in their 
respective topological sectors. We will call these states the RK states. 
This is the first example in these lectures of a Resonating Valence 
Bond wave-function. 

• It has been shown by Kohmoto and Shapir ^H2], that the spin-spin 
correlations in this state decrease exponentially. 

• An important property: any dimer correlation functions in the RK 
state can be computed from an exact mapping to the classical statis- 
tical problem of dimer coverings first solved by M. E. Fisher and J. 
Stephenson jH^]. From this work one can conclude that the dimer- 
dimer correlation functions at y = 1, decreases algebraically with dis- 
tance (as r~^). This property implies that the first excitations above 
the ground-states are gapless. 

• On the basis of the continuity in the energy between the staggered 
phase and the RK states, one may speculate that the RK point is the 
quantum critical end of the staggered VBC phase. But the excitations 
of the staggered VBC are non local and have energy of order 0{N^^'^) 
in the -j — > oo limit. To sustain the above point of view one should 
explain how the kinetic term can dress these excitations so that they 
become gapless when -j — > 1. In fact the more probable hypothesis is 
a first order phase transition between the RK phase and the staggered 
one. Such a question could perhaps be answered with Monte Carlo 
simulations. 

• The ground-state wave-function at this RK point has a property which 
is considered as constitutive of a RVB spin liquid: that is resonances 
between all dimer coverings. It must be underlined here that this res- 
onance phenomenon at the RK point does not bring any stabilization 
of the equal amplitude superposition ground-state when compared to 
the neighboring staggered VBC phase. 

• RK then argue that for j- < 1, there is, separated by a first order phase 
transition, a new phase which might be a "true" resonating Valence 
Bond Spin Liquid ^. The characterization of this phase is for the 
moment rather loose: RK argument is variational and rather week. 
The first calculation by Sachdev on a 36 lattice i84j, gives evidence 

^i.e. a phase where the resonances between different dimer coverings are essential to 
its stabilization and are so important that there is a gap to the first excitations and any 
correlation functions: either spin-spin, dimer-dimer or higher order plaquette-plaquette 
correlation functions decrease exponentially with distance. 
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for a VBC columnar state for -j < 0.5 and not a clear conclusion 
nearer from the RK point. Extending the calculations to 64 sites, 
and using various estimators, Leung and co-workers J85] estimated 
that long range columnar order probably exists up to the RK point, 
with the restriction that up to -j ~ —0.2 the order is very plausibly 
purely columnar, whereas in the range —0.2 < ^ < 1 the order could 
reduce to a 4-spin S=0 plaquette order ^. It seems nevertheless widely 
admitted |86[ 18 7j that this model has crystalline order everywhere 
except at the RK critical point. 

In view of these results for the QHCD model on the square lattice, of 
most studies on the Ji — J2 SU{2) model, and of the SU (N) studies on 
the same lattice [HHl En]; one may be tempted to conclude that VBC is 
the paradigm of the quantum ground-state on square and possibly bipartite 
lattices. This might be an escapable assumption |511l20j . but the fact is that 
the triangular based lattices (see next chapters) seem much more favorable 
to Resonating Valence Bond Spin Liquids. 



*This conclusion is not consistent with the degeneracy the authors claim for the ground- 
state. The nature of the phase for -j < 1 remains an interesting and open question: 
interesting but technically difficult. The same kind of difficulty is present in the study 
of the Ji — J2 model on the square lattice for J2/J1 ~ 0.5. At this point of maximum 
frustration, Neel order is destroyed but the exact nature of the phase is uncertain: colum- 
nar order |46l 148) . 4-spin plaquette order [471 15()| or RVB spin liquid [51|? In view of 
exact spectra for sizes up to N=36, it seems that the 4-spin S—O plaquette order is the 
less plausible (because the k — (tt, tt) states necessary for the 2-fold symmetry breaking 
of this state is very high in the spectrum). We expect a 4-fold symmetry breaking in 
the columnar state as well as in the RVB state and refs. therein). The gaps from 

the ground-state to the plausible candidates for these 4-fold symmetry breakings are still 
very large in the N=36 sample. We are thus lead to conclude that the N=36 sample is 
too small to give informative issue on the dilemma: columnar state or RVB state. This 
strongly weakens the variational argument of ref. |51| . 
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Chapter 5 



Resonating Valence Bond 
Spin Liquid (Type I) 



The Resonating Valence Bond Spin Liquid is a quantum concept introduced 
in 1973 by P. W. Anderson [T3], following the line of thought of Linus Paul- 
ing for molecules. When the semi-classical Neel states or simple dimer 
covering solutions are very far to satisfy each individual bond, Anderson 
speculated that the macroscopic system could take advantage of the quan- 
tum resonances between the exponential number of dimer coverings to lower 
its ground-state energy. Such states have no long range order whence the 
name of Spin Liquid, quantum resonances between the exponential number 
of equivalent dimer coverings are essential: it is a Resonating Valence Bond 
Spin Liquid (abbreviated as RVB Spin Liquid or RVBSL in the following). 

5.1 Introduction: short range versus long range 
Resonating Valence Bond wave- functions 

Resonating Valence Bond wave-functions encompass a large class of wave- 
functions beyond the equal amplitude superposition of next neighbor dimer 
coverings that we encounter at the RK point in the last chapter. 

It is easy to verify that the whole set of dimer coverings (without any 
restriction on the length of the dimers) is an overcomplete basis of the S = 
subspace of the spin system (compare the numbers of these coverings to the 
size of the S = subspace for a N site lattice) ^. 

Let us suppose that we have designed a family £ of linearly independent 
dimer coverings Cj, a general RVB wave- function will be written as: 



For large enough sizes the next neighbor coverings form a hnearly independent family. 





62 



where \Ci > are products of dimer wave- functions (with a sign conventionahy 
fixed, respecting the lattice topology). 

In variational calculations, one generally use restricted forms of Eq. 15.11 
where the amplitude A{Ci) of a given configuration Ci is written as the 
product of amplitudes h{k,l) for each dimer {k,l) present in Cj. 

Two situations have been studied: 

i) either long range RVB wave functions where the function h{k, I) de- 
pends algebraically on the distance r^i between sites k and I (at least for 
large distances): 

hik,l) = — (5.2) 

^kl 

Liang, Dougot et Anderson have shown that such wave functions have 
Neel long range order in the Heisenberg model on the square lattice if a < 5 
and no Neel long range order for o" > 5 ^. Capriotti and co-workers |51j 
have used a p-wave BCS wave-function for the Ji — J2 model on the square 
lattice, which has no long range order in dimers. 

i) or the short range Valence Bond w.-f. where the amplitudes h{k, I) 
are not necessarily strictly restricted to next neighbors but decrease at least 
exponentially with distance (most of the following is concerned with that 
kind of wave functions). By construction such functions cannot describe 
Neel long range order, as Neel order has long range correlations between 
spins on the same sublattices. As we have seen in the previous chapter, this 
family encompasses the quantum critical behavior of the QHCD model on 
the square lattice. We can equally describe in this basis the Valence Bond 
Crystals, which are characterized by dominant amplitudes associated to the 
simple symmetry breaking configurations. Many properties of these wave- 
functions have been studied theoretically ( jHD ESI ESI EOl 1^ and references 
therein), we will see some of them in the following. 

In this chapter we will first describe with some length the properties 
of the QHCD model on the triangular lattice, to compare to the solution 
of the same model on the square lattice. We will then move to the Multi- 
Spin Exchange Hamiltonian on the same lattice, which is the first SU{2) 
model exhibiting a "true" resonating Valence Bond Spin Liquid. A special 
attention will be given to the topological degeneracy of the ground-state, 
and to the existence of deconfined spin- 1/2 "spinous" excitations, which is 
the most important experimental signature of a RVB Spin Liquid state. We 
will close the chapter by a small bibliography on gauge theory approaches 
that have been dealing with the same physical problem. 

^They equally show that the difference in energy of those different wave- functions are 
extremely tiny 
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5.2 The Quantum Hard Core Dimer model on the 
triangular lattice 

The QHCD model on the triangular lattice has been studied by Moessner 
and Sondhi in 2001 jST], when they realized that the dimer-dimer corre- 
lation function on this lattice was not algebraically decreasing as on the 
square lattice but exponentially decreasing with distance. The model on the 
triangular lattice comprises the same ingredients as on the square lattice: a 
potential energy term between parallel pairs of dimers and a kinetic energy 
term which does a cyclic permutation of parallel dimers on 4-spin plaquettes 
(involving two triangular units). 



Plaquettes 



+ 



(5.3) 

The sum over plaquettes runs on the three kinds of plaquettes with orien- 
tations at 60 degrees from each other. 

J can be assumed to be positive ^: it enforces resonance effects, V can be 
positive (repulsion between dimers) or attractive. The conditions of validity 
are the same as those of the model on the square lattice: it is supposed 
that the spin gap is large enough so that the first excitations are in the 
singlet sector. Insofar as the spin gap is large, the spin-spin correlations 
are short range which is consistent with the restriction to the subspace of 
nearest-neighbor Valence Bonds. 

The properties of the lattice affect the properties of the QHCD model 
on two central points: 

• In the triangular case due to the higher entanglement of the lattice 
with the two-dimer terms, there is only 4 different topological sectors 
classified according to the parity of the winding numbers: (even, even), 
(even, odd), (odd, even), (odd, odd). (The dimer-flip term can change 
the winding numbers, not their parities). 

• At temperature much larger than J and V, the square lattice problem 
has algebraically decreasing dimer-dimer correlations, whereas on the 
triangular lattice these correlations decrease exponentially |87j . 

As in the square lattice case, at the point J = V the model is exactly 
solvable. The ground-states here too are the equal amplitude superpositions 
of all dimer coverings in each topological sector. 

Demons: A lower bond for the ground-state energy is given by a min- 
imization of the Hamiltonian on each plaquette individually. If the given 
plaquette has no parallel dimer (non flippable plaquette), its energy is zero 

^This was not obvious a priori and is important in the following as it insures that the 
Hamiltonian is positive semi-definite for V > J > 
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and if it has parallel dimers its potential energy is V and its kinetic en- 
ergy is < — J. We can thus write a lower bond for the energy of the global 
system, which is proportional to the number of flippable plaquettes njnp, 
as min [0, {V — J)nfiip]. At V = J the equal amplitude wave-functions in 
each topological sub-sector saturate this lower bond. As TCdimer is positive 
semi-definite at this point, the equal amplitude wave-functions are thus the 
unique ground-states of the problem. 

Contrary to the case of the square lattice, the degeneracy of this RVB 
subspace is only 4 on the triangular lattice (whereas it is of order N in the 
square lattice case). As in the square lattice case, these RVB states are 
degenerate with the 6 staggered configurations, which are ground-states for 
any V/J > 1 (see Fig. lOT) . 

The sum over all configurations of the equal amplitude wave-functions is 
equivalent to the classical dimer problem (up to the question of the staggered 
phase which has a negligible statistical weight in the problem): thus the 
dimer-dimer correlations decrease exponentially with distance at the point 
J = V. It's the description of a true RVB Spin Liquid phase with exponen- 
tially decreasing spin-spin and dimer-dimer correlations (we thus expect a 
gap in the singlet sector), and translational invariance of the ground-state 
(all four topological ground-states, with equal amplitude wave- functions are 
in the same k = (0,0) sector of the impulsion |81j). 

Monte-Carlo simulations |SZ| have shown that this phase extents at least 
in the range 2/3 < V/ J < 1. It terminates at V/J = 1, with a first order 
transition to the staggered phase (seen in the Monte Carlo simulations as 
hysteretic behaviors). The dimer-dimer correlation function is very short 
range in all the above-mentioned range of parameter, and very weakly de- 
pendent on temperature, which is suggestive of a gap in the spectrum. 

We will not comment on the rest of the phase diagram, as it is not 
relevant to our main point here. It is described in Fig. 15.11 (For more 
details, see the original paper [H7]). 

Spinons: As noted above the RVB phase has a gap to collective exci- 
tations, which is equally true of VBC. The major difference insofar between 
VBC and RVB Spin Liquids is the existence in this new quantum phase of 
deconfined spin-1/2 excitations: the spinons. If you break a Valence Bond in 
a VBC phase and try to separate the two single spins from each other the en- 
ergy of the system increases as the length of the string of misaligned dimers 
which appears between the two single spins (take as an example the stag- 
gered or the columnar phase of the QHCD). This creates an elastic restoring 
force which binds the two spin-1/2 together: in such a Valence Bond Crys- 
tals excitations have always an integer spin (AS = or 1). We suspect that 
in the RVB Spin Liquid state, where the correlations between local oper- 
ators are short range and any disordered configuration as probable as an 
other, the restoring force between two single spins beyond a certain distance 
will be negligible and the spin-1/2 ("spinons") will be deconfined. A simple 
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Figure 5.1: The phase diagram of the Quantum Hard Core Dimer problem 
on the triangular lattice 



verification can be done on the equal amplitude states of the Quantum Hard 
Core Dimer model whatever the lattice: spinous do not interact beyond one 
lattice step. One expects this property to extend in all the RVBSL phase 
at T=0. At high temperature the classical square lattice is known to be 
logarithmically confining j^. Moessner and Sondhi have checked that the 
triangular lattice is not confining. 

The existence of deconfined spin-1/2 excitations, and thus of a continuum 
of excitations just above the gap is the main experimental signature of the 
2-dimensional Valence Bond Spin Liquid. It was recently claimed that this 
continuum of excitations has been observed in CS2CUCI4 which is supposed 
to be a two dimensional magnet |94j . 

Two more complex evidences of RVB Spin Liquids phases had been ob- 
tained before the discovery of this simple toy model: the first in 1992 in 
a large N, Sp{N) analysis of Sachdev j25J) the second in a more realistic 
SU{2) spin model by Misguich and coworkers [961 197j . This last work will 
be the object of the next section. 



5.3 The MSE model or Ring Exchange model on 
the triangular lattice 

The multiple-spin exchange model (called MSE in the following) was first 
introduced by Thouless [HHl to describe the nuclear magnetism of three- 
dimensional solid He^ and by Herring 100 for the Wigner crystal. It 
is an effective Hamiltonian which governs the spin degrees of freedom in a 
crystal of fermions. The Hamiltonian is a sum of permutations which ex- 
change the spin variables along rings of neighboring sites. It is now largely 
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believed that MSE interactions on the triangular lattice also describe the 
magnetism of solid He^ mono-layers adsorbed on graphite jlUlll96MlU2j and 
that it could be a good description of the two dimensional Wigner crys- 
tal of electrons jlU3j . In He^, exchange terms including up to 6 spins are 
present |inij . Recent discussions equally concern the strength and impor- 
tance of the 4-spin exchange term in La2Cu04 jl04t ?, ?,|S2]- 

Here we will only focus on 2- and 4-spin interactions which constitute 
the minimal MSE model where a short-range RVB ground-state is predicted 
from exact diagonalizations ^7\. The Hamiltonian, which is also called by 
some authors the Ring Exchange model, reads: 

The first sum runs over all pairs of nearest neighbors on the triangular 
lattice and Pij exchanges the spins between the two sites i and j. The second 
sum runs over all the 4-sites plaquettes and Pijki is a cyclic permutation 
around the plaquette. The 2-spin exchange is equivalent to the Heisenberg 
interaction since Pij = 2Si ■ Sj + 1/2, but the four-spin term contains terms 
involving 2 and 4 spins and makes the model a highly frustrated one. 

The general phase diagram of this model is given in Fig. 15.21 

We will now focus on the phase described as "Spin Liquid I" in Fig. 15.21 
and more precisely on the point J2 < and J2/ Ja — —2 which has been stud- 
ied extensively by means of exact diagonalizations up to = 36 sites (HZ] 
(this is a good qualitative description of the low-density solid He^ films). 

Finite size effects on the spin gap and energy per spin are displayed in 
Fig. 01 

These data point to a spin-gapped phase with a short correlation length 
(of the order of a few lattice steps) and a spin gap of order 1. 
Three properties should be emphasized: 

• No sign of a VBC could be found. All correlations functions spin-spin, 
dimer-dimer, 4- and 6-spin plaquette-plaquette seem short range j97| 
llU8j and consequently all susceptibilities associated to local observ- 
ables are zero in the ground-state of the MSE spin liquid (see ref . [HH ) • 

• The ground-state displays at the thermodynamic limit a degeneracy 
that has been shown to be purely of topological origin j81j . 

• The system probably supports unconfined spinous ^. 

*On the basis of too small samples we had concluded in our 1999 paper that spinous 
were probably confined (which was a bit unpleasant and contradictory with the existence 
of a topological degeneracy and the absence of any order in a local order parameter.) 
Extending the calculations to larger sizes up to A'' = 33, we have now results that clearly 
point to deconfined spinous for large enough distances. 



67 



Class 

ft 


ical Phas 


>e Diagram 


Kuboe 

■ 


tMomoi 97,98 


-4 


uuud 
I +... 




1. ~ 


Neel 120° 

5.8 

■ 



QU£ 

ft 


intum Phase Diagram 

(^^Spin Liquid!^ 


.iquid II 

lisguich '00 


t 




5 Misguich 98,99 +1 


Spin i 

LiMing, A 
+ 


jf/'j, 

10 



Figure 5.2: The phase diagram of the ring exchange model fEa. I5.4j) on the 
triangular lattice [Tn7|l97]. 
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Figure 5.3: Finite size effects on the spin gap (top graph) and energy per 
spin (bottom graph) in the Ring Exchange model (Eq. 15. 4j) for J2 = —2, 
J4 = 1. Ref. inn and unpubhshed results. Samples with an odd number of 
sites are indicated by crosses. The full squares are for even samples with 
the full symmetry of the infinite lattice and open squares for even samples 
with lower spatial symmetries. Study of the energy per site (bottom figure) 
is specially interesting. Shape effects are still important for sizes as large 
as 32 but is appears clearly that the energy of the most frustrated samples 
(odd number of sites: crosses) is converging to the same limit as the energy 
of the unfrustrated ones. This is a good indication that the larger sizes that 
we have considered should allow significant qualitative conclusions. 
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5.3.1 Topological degeneracy: 

This subject is fully developed in ref. and illustrated in Fig. 15.41 In 
this figure the first graph (top left) shows the low lying singlet states for 
the N=36 sample. On this graph one sees that these singlet levels appear 
as multiplets (the black symbol is one time degenerate and the red symbol 
has a 3- fold degeneracy). With the system size the "red" levels collapse to 
the black ones exponentially fast, with a characteristic length which is of 
the order of 0.6 lattice step (finite size scaling in the bottom right figure): 
this is the degeneracy, that we argue to be of topological origin. 

In short the topological degeneracy can be understood using two argu- 
ments: i) the wave- functions describing the MSE spin liquid can be classified 
as short range RVB wave functions (all the correlations functions in local 
observables are short range), ii) the Hamiltonian is a local operator insen- 
sitive to a global property as the parity of the winding number. In other 
words for large enough sizes, it is possible to locally optimize the energy and 
the result does not depend on the topological sector where it is done. As 
there is 4 topological sectors for a triangular lattice on a two torus, we have 
there the origin of the 4-fold degeneracy. 

On lattices with an odd number of rows it is possible to transform one 
topological sector in another by a 2tt twist of the boundary conditions along 
the even direction (Fig 15.41 top-right graph). ( This is equivalent to the 
introduction of a quantum of fictitious flux through the torus: detailed 
explanations can be found in ref. j81j^. In such an operation the global 
spectrum of the Hamiltonian is unchanged but eigen-states are not mapped 
on themselves: the momentum of the states is translated by (vr, 0) in such 
an operation. Oshikawa jl09j concluded that such systems should exhibit a 
doubling of the Brillouin zone and the spontaneous symmetry breaking of 
a Valence Bond Crystal. This is an incorrect speculation as we will explain 
in section 5.5. A simple counter-example is given by the HCQD model: the 
ground-state in the different topological subspaces for (even, even) samples 
are all in the k = (0, 0) sector of the momentum. There is thus an alternation 
of the quantum numbers of the degenerate ground-state multiplicity when 
going from (even, even) to (even, odd) samples whatever the system size: 
this indeed is inconsistent with a VBC spontaneous symmetry breaking in 
the thermodynamic limit. In the same line, in the MSE model, the evolution 
with the system size of the quantum numbers associated to the point group 
equally shows that this topological degeneracy is by no way associated to 
some symmetry breaking of the lattice (81) . 

A last comment relative to the possibility of using this degeneracy to 
produce quantum bits protected from decoherence effects 
At first sight the idea is attractive: due to the absence of long range order 
in any local variable the local susceptibilities vanish in the thermodynamic 
limit (see ref. ^T\) and one expects such quantum bits to be insensitive to 
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Figure 5.4: Topological degeneracy in the Ring Exchange model (Eq. 15. 4j) 
for J2 = —2, J4 = 1. Ref. [HZj and unpublished results. The figure at the top 
left shows the low lying levels in the singlet and triplet sectors for the sample 
of size N = 36. Notice that the singlet gap is already much smaller than the 
triplet gap. The figure at the bottom left gives the finite size scaling of the 
spin gap. The figure at the bottom right shows how the gap in the singlet 
sector is closing as a function of the size. It is exponentially decreasing with 
the linear size of the lattice, and the correlation length is about 0.6 lattice 
step. Such a law correctly describes the finite size scaling of the spin gap 
of unfrustrated samples (bottom left graph), with an estimate of the spin 
gap ~ 0.8 lb 0.1. The figure in the top right shows the effect of an adiabatic 
twist of the boundary conditions in exchanging the two quasi-degenerate 
topological levels. 
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Figure 5.5: Spinons binding energy in the Ring Exchange model (Eq. I5.4j) 
for J2 = —2, J4 = 1. One sample only has been used for the sizes 28, 29, 
30; this precludes an estimation of the uncertainty due to the sample form 
for intermediate results. Notwithstanding this, the displayed results do not 
infirm the conclusion of a zero energy binding for large enough sizes. The 
red curve is an imperfect and arbitrary polynomial fit to indicate a general 
tendency (unpublished results, work in progress). 

any local cause of decoherence. But the charm of this property is to be paid 
by highly non trivial, if not impossible, writing and reading of the state of 
the quantum bit (which would imply manipulation of gauge fields ..). 

5.3.2 Unconfined spinons? 

To estimate the confinement energy of two separated single spins (see Fig ]5.5() , 
we use the following arguments: 

ij A sample of 2N spins has a ground-state energy which is essentially: 

E{2N,S = f)) = 2Ne^ + 0{^), (5.5) 

v N 

the correction term to the N ^ 00 limit is at most 0{-^) if there is long 

range order in spin-spin correlations and is more plausibly 0{ ^^^^^^^^ ) in 
the present case. 

ii) A sample withe 2N + 1 sites can accommodate N valence bonds -|- 
one single spin, we can write its ground-state energy as: 

E{2N + 1, 5 = 1/2) = (27V + 1) eoo + ^sp^non + ■ ■ •) (5.6) 

V N 

where Aspinon measures the energy gap for the creation of one spinon. 

in) If we now look to the first AS" = 1 excitation of an even sample, 
we expect (in the hypothesis where spinons are the first excitations of this 
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system) that: 



E{2N, 1st exc. lev.) = 2Neoo + 2Aspinon + Ebinding + • • •) (5-7) 

where Ehinding is the binding energy of two spinons. If Eififiding is negative in 
the thermodynamic hmit then the spinons wih be confined, the first excita- 
tions will be integer spin excitations as expected in a Valence Bond Cristal. 
If Ebinding with Aspinon — > / with increasing size, we then ex- 
pect the spinons to be unconfined and the first excitations of the model are 
fractionalized [TT^ ITT:?! [TT?1 ITT^ [TTTH ITTTl ITTRj 

Using these three equations fEas. 15.51 15.61 I^TTj) for consecutive sizes one 
obtains an estimate of the binding energy for each sample size. The results 
are shown in Fig. 15.51 They give a positive indication in favor of unbound 
spinons in the MSE model. 

5.4 RVB Spin Liquids in other spin models 

We suspect that RVB Spin Liquids could be observed in other spin models. 
As already discussed in the previous chapter, the case of the Ji — J2 model 
on the square lattice is still debated, the correlation length being probably 
larger than the largest sizes actually available. The Ji — J2 model on the 
hexagonal lattice may have spin liquid phases for J2/ Ji ~ 0.3 and around the 
point with ferromagnetic J2 and J2/J1 = —0.25 |20j . Here too, and contrary 
to the MSE model the range of parameters where a spin liquid phase might 
appear is relatively small, the local S=0 objects probably extend over a few 
lattice cells, and as a consequence the gaps are rather small and the shape 
effects a bit chaotic for the available sizes. 

An explanation of the robustness of the short-range RVB phase in the 
MSE model can be guessed from the analogy between multiple-spin interac- 
tions and QHCD models. From the analysis of QHCD models we understand 
that RVB phases are possible when VBC are energetically unstable. Colum- 
nar VBC are stabilized by strong parallel dimer attraction and staggered |HZj 
VBC appear when the repulsion between these parallel dimers is strong. In 
between, an RVB phase can arise ^. From this point of view, increasing the 
role of the kinetic term of the model is essential: it is exactly the role of the 
four-spin ring exchange term of the MSE model (Eg. 15.4)) 108]. 

The role of the triangular lattice (or of a preferred triangular sublattice 
as in the case of the Ji — J2 model on the hexagonal lattice) should probably 
also be emphasized: J.C. Domenge's preliminary work on the MSE model on 
the square lattice points to a smaller extent of the Spin Liquid phases j32j . 

^In Ref. [1181 an RVB state is selected by introducing defects in the lattice in order to 
destabilize the competing VBC states. 
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Figure 5.6: 2-torus with one cut A. 

5.5 Short range RVB: topological degeneracy and 
absence of symmetry breaking 

In this section we develop two important properties of type I short range 
RVB Spin Liquids, in a more general point of view than in the preceding 
sections and independently of a specific Hamiltonian, with the only restric- 
tion that it is a short range operator and that the ground-state and first 
excited levels can be described with resonant Valence Bond superpositions 
of short range dimers (not exclusively first neighbors). 

(For simplicity this section includes some parts of the paper by Misguich 
et al. [Hl\ but some demonstrations and examples, which should not alter the 
general understanding are omitted and in some cases the reader is referred 
to the original work to complete the picture.) 

5.5.1 Topological degeneracy of the ground-state multiplic- 
ity of a type I RVB with half integer spin in the unit 
cell 

Definition of the topological sectors 

Let us draw a cut A encircling the torus created by periodic boundary con- 
ditions (see Fig. 15.6)1 . This hyper-surface of dimension d — 1 cuts bonds of 
the lattice but there is no site sitting on it. The position of the cut is ar- 
bitrary. The family of nearest-neighbor dimer coverings can be decomposed 
into two subspaces £^ depending on the parity Ha of the number of dimers 
crossing the cut A ^. By considering a set of d cuts Aj=i ...^^ encircling the 
torus in all possible directions one obtains 2*^ families of dimer covering. 

Any movement of dimers can be represented as a set of closed loops 
around which dimers are shifted in a cyclic way. A local operator will only 
generate contractible loops which will cross each cut a even number of times. 
The number of dimers crossing the cut can therefore only be changed by an 
even integer and the parities are unchanged. 

^The definition of Rokhsar and Kivelson |79| that we used before is equivalent to this 
one but less practical in the present context. 
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This property remains true as long as one works in a subspace where 
the dimer lengths are smaller than the linear system size, that is when the 
topological sector are well defined (if a dimer length is half the linear of the 
system one cannot decide by which side it goes). On the other hand, we 
checked on the triangular (resp. Kagome) lattice that these 4 sectors are 
the only topological sectors: local 3- (resp 4-) dimer moves can be used to 
transform a configuration of a given sector into any other configuration of 
that sector. 

These subspaces are orthogonal in the thermodynamic limit jll9| . The 
graph of the scalar product < c+|c^ > of two dimer configurations belonging 
to different subspaces E'^ and E~ has at least one long loop encircling the 
torus in the direction. When goes to infinity this contribution to 
the scalar product is smaller than 2~^^/^. Consider two normalized vectors 
1^'"'') and |^'~) belonging to two different sectors: 



Because of the exponential number of dimer coverings in each subspace it is 
not obvious that 1^'"'') and l^*") are orthogonal in the thermodynamic limit. 
But it is nevertheless the case [5^, and: 



In the following, unless explicitly mentioned, we consider the 2D case 
for simplicity but most of the topological arguments about dimer covering 
immediately extend to higher dimensions. 

Two- fold degeneracy in even x odd samples 

In the special case of tori with an odd number of rows (and an odd number 
of spin- 1 per crystallographic unit cell), one step translation along the x axis 
(called 7^ in the following) maps E'^ on E~ and reversely. Some point-group 
symmetry can also do this job. A vr rotation about a lattice site nearby the 
cut (called TLt^ in the following) has the same effect. If the cut is chosen 
parallel to a symmetry axis of the sample, a refiection with respect to this 
axis (called Yjy in the following) will equally map E~^ on E~ and reversely. 

All these symmetry operations isolate a single column C of lattice sites 
between A and its transform A'. In that case columns have an odd number 
of sites and an odd number of dimers must connect some sites inside C with 
sites outside C. Therefore IIa differs from 11^' and the two subspaces E^ 
and E~ are exchanged. 

For a large enough system these two sectors are 1) orthogonal, 2) uncou- 
pled by any local Hamiltonian and 3) exchanged by symmetry operations 
(even X odd). This is enough to insures that they have the same spectrum. 






(5.9) 



75 



irrespectively of the physics of the model, provided it can be described in 
the short-range dimer space. In fact, quantum numbers of these doublets of 
degenerate states are fixed by symmetry. 

We decompose an eigenstate \^po > on the two topological subspaces 
defined relatively to the cut A: 

lV'o) = |^o^> + |V'o") (5.10) 

where |V'o^) belong respectively to the sets \tpo), as an eigenstate of the 
Hamiltonian with periodic boundary conditions, belongs to an irreducible 
representation of the translation group. In the following we will also assume 
an TZt^ and invariance of the Hamiltonian and that, for simplicity, IV'o > 
transforms under a one-dimensional representation under TZj^ and T,y. 



'r.|V'o) = e^'^°-"|^o) 

S^lV'o)=^o'lV'o) (5.11) 

In the thermodynamic limit \4'o) and \^o) orthogonal and 7^., TZ^^ and 
map on £~ and reversely: 

Let us now build the variational state: 



|V'i,a) = |^o^>- IV'o"). (5.13) 

Eqs. (f^TT^ imply: 

% I^Pia) = IV'i.a) 
T^n |V'i,a) = -Po IV'la) 

S^|V'i,a) = -'^oIV'i,a) (5.14) 

\iPi.a) has thus a wave- vector ki, a rotation quantum number and a 
reflection quantum number af related to the quantum numbers of {iPq) by 
the relations: 

ki = ko + (vr, 0) 
PJ = -Po 

al = -al (5.15) 

It is thus a state orthogonal to the ground-state (even on a finite-size system 
where the topological sectors are not rigorously orthogonal). 
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Since any local Hamiltonian has exponentially vanishing matrix elements 
between different sectors we have 



and IV'i.A > is thus degenerate with the absolute ground-state, their sym- 
metries being related by relations (|5.15)) . 

Degeneracies in even x even samples 

As remarked by Bonesteel |119j on a square lattice a 7r/2 rotation exchanges 
sector — ) and sector (— , -|-) but sectors (— , — ) and (-I-, -|-) remain inequiv- 
alent. A similar phenomenon occurs on the triangular lattice where 27r/3 
rotations permute cyclically 3 of the 4 sectors. As a result there are equally 
in even x even samples exact degeneracies due to these mappings SI , . Read- 
ers interested in the details of these degeneracies are referred to the original 
paper [HT] . 

4-fold degeneracy in RVBSL phases 

The numerical data on the type I Spin Liquid phase of the MSE and QHCD 
models jH^ suggest that the ground-state degeneracy is 4 in such a phase 
whatever may be the shape of the sample. 

There is no global mathematical proof of this property but the following 
physical arguments make it extremely plausible. 

Let us make the following assumptions: 

a) The ground-state can be described in a short-ranged dimer basis, b) 
All n-dimer correlations (n = 2, 3, • • •) are short-range and the corresponding 
correlation lengths are bounded, c) The Hamiltonian is local. 

From hypotheses a) and c) it is clear that for a large enough system the 
four topological sectors are not mixed in the ground-state and the spectrum 
can be computed separately in each sector. As seen before, we do not have 
any symmetry operation which connects all four sectors (these operations 
connect only two or three of the topological subspaces, depending on the 
geometry of the sample) and we need a physical argument to explain why 
energies should be the same in each sector (in the thermodynamic limit)? 

Because of their topological nature, it is not possible to determine to 
which sector a dimer configuration belongs by looking only at a finite area 
of the system. In other words, any dimer configuration defined over a large 
but finite part of the system can be equally realized in all sectors. The 
Hilbert space available to the system is the same over any finite region of 
the system. In the absence of any form of long-range order the system can 
therefore optimize all its correlations with an arbitrary high accuracy equally 
well in each sector. At this point we can only conclude that the four sector 




(5.16) 
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will have the same energy density and we cannot exclude the existence of a 
gap between the different topological sectors. 

However the numerical results obtained in the QHCD and MSE models 
indicate that it is not the case and that the four ground-states have asymp- 
totically the same total energy. We think that this should be true for a 
general short-range RVBSL. 

Miscellaneous remarks on the RVB ground-state degeneracy 

— Dimers and twist operator. The variational state |V'i,a) can be deduced 
simply from |^o) by changing the sign of the dimers crossing the cut A. 
Such an operation can also been seen as a 2tt twist of the spins of column 0. 
From the physical point of view the reason why IV'i.a) has asymptotically 
the same energy as becomes clear: in the absence of stiffness (~ ab- 
sence of sensitivity to a boundary twist in the thermodynamic limit) and of 
long-range spin-spin correlations, the perturbation induced by the bound- 
ary condition cannot propagate and does not change the energy of the initial 
state: its only effect is to change the relative phases of the different topo- 
logical components of the wave function, and consequently the momentum 
and space symmetry quantum numbers of the initial state of the even-odd 
samples. 

— Fractionalization and topological degeneracy. To our knowledge all 
present theoretical descriptions of fractionalized excitations in 2D magnets 
or related problems jl2Ul 11151 IT^ l . 112, 118 (we should also mention topo- 
logical properties of Laughlin's wave function for fractional quantum Hall 
effect |1221 1123] ) imply topological ground-state degeneracies. In such pic- 
tures, the physical operation which transforms a ground-state into another 
is the virtual creation of a pair of spinous (by dimer breaking) followed by its 
annihilation after the circulation of one of them around the torus. In such a 
process a vr phase-shift is introduced between the topological sectors (as in 
the above recipe). For samples with an odd number of rows this operation 
connects eigenstates with different k vectors (and space quantum numbers) 
as described in Eqs. ()5.15|) . 

— Numerical studies with a different topology. An interesting check of 
the pure topological nature of this degeneracy could be obtained by studying 
the problem no more on a torus but on a surface with a different genus. On 
a sphere we expect an absence of degeneracy. Unfortunately if a lattice 
can be represented on an infinite plane, both the number of links L and 
plaquettes P depend linearly on the number of sites and Euler's relation 
P — L-\-N = 2 — G constrains the genius G to be 2 ! The torus is the only 
possible topology if we require a full translation invariance in both directions. 
In a recent work, loffe et al. |lllj have studied the absence of sensitivity to 
disorder as an evidence for topological phenomenon in the liquid phase of 
the QHCD model on the triangular lattice. They also used open boundary 
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conditions to modify the topology of the system and argue in that case that 
the low-energy spectrum is free of edge states which could hide the actual 
ground-state degeneracy. 

— Example of RVB phase with 2 spins in the unit cell. A spin liquid 
state, seeming very similar to the state observed in the MSE model on the 
triangular lattice, has been observed in the Ji — J2 model on the hexagonal 
lattice [20] for Ji = — 1, J2 = 0.3. No quasi-degeneracy of the ground-state 
has been noticed. It should be remarked that in this system there are 2 
spins I per crystallographic unit cell and no degeneracy is expected on the 
basis of the topological arguments. 

5.5.2 Symmetry breaking in gapped phases 

From the mathematical point of view, ground-state wave functions that 
break one-step translations or space group symmetries can be built from 
linear combinations of the degenerate ground-states of the even- odd sam- 
ples. In a completely equivalent way, ground-state wave functions that 
break rotation symmetry can be built in even x even samples. One could 
thus superficially conclude that spontaneous symmetry breaking is possible 
in RVBSL, we will show below that this assumption is false. 

There are many features which show that this degeneracy property is a 
subtle one, both from the mathematical and physical viewpoints. 

The possible alternation of the spatial properties of the low-lying excita- 
tions with the parity of the number of rows of the sample (as observed in the 
QHCD model on the triangular lattice) is a first difficulty. The degeneracy 
of the RVBSL is in fact quite different from that appearing in a VBC. We 
do not expect the VBC ground-state degeneracy to depend on the genus of 
the sample, as the RVBSL does. 

Prom the physical point of view also the two situations are quite differ- 
ent. An infinitesimal symmetry breaking perturbation is able to select one 
symmetry breaking ground-state of the VBC, but as we will show below this 
is impossible in the RVBSL. 

Let us call A the extensive non-diagonal observable appearing in the 
VBC in the thermodynamic limit. On a columnar VBC modulated in the 
u direction, this observable is: 

N 

A = Y,e''^'"''Ps=o{rj,rj + u) (5.17) 
i=i 

where Ps=o(rj, rj + u) is the projector on the singlet state of two neighboring 
spins. A connects eigen-states with wave-vector ko to states with wave- 
vector ko + Ki . 

On a finite size sample, with periodic boundary conditions, the expecta- 
tion value of A is zero in any eigenstate, but < A'^ > could be non zero. If 
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the order parameter V defined by: 

=<^l:g.s.\A^A\il:g.s.> /N^ (5.18) 

does not vanish in the thermodynamic hmit, the system has columnar dimer 
long range order with wave vector Ki. 

Let us now consider a perturbation of the Hamiltonian: 

Hs = Ho- {6A + h.c.) . (5.19) 

At T=0, the intensive linear response on the observable A is measured by 
the susceptibility: 

X = ^ < i^g.s.lA^ rr \ -^li^a.s. > (5.20) 
This susceptibility is bounded from below P^lj : 



/ 

where / is the oscillator strength: 



< X (5.21) 



f = ^{i;g.s.\[A,[Ho,A]]\i;g.s.) (5.22) 



The demonstration uses the properties of the spectral decomposition asso- 
ciated to the operator A: 



S{Lo) = ^Y.\<i>g.s.\A\n>f6{io-LOn) (5.23) 
where u;„ = En- Eg,, 



P2 = 1 1 S{u)du; (5.24) 
Using the Cauchy Schwartz inequality one obtains: 

<^ j ujS{uj)duj j uj-^S{uj)duo (5.25) 

j ujS{uj)duj = f/2 (5.26) 



where 



uj-'S{u)duj = x/2 (5.27) 

which proves inequality ()5.21|) . For a short range Hamiltonian the oscillator 
strength / is 0{1) and inequality (|5.21|) implies that the T=0 susceptibility 
associated to a finite order parameter diverges at least as the square of the 
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sample size: any infinitesimal symmetry breaking perturbation will select a 
symmetry breaking state. 

We will now show that for a RVBSL, where all the correlations functions 
are short-ranged with correlation lengths bounded by the susceptibilities 
of the medium remain finite in the thermodynamic limit. To do so we 
distinguish in Eq. l5.2()l the contributions from the quasi-degenerate states of 
the topological multiplet (called \ai >) from the contribution of the other 
states of the spectrum, above the physical gap A. We thus obtain the 
following upper bound for the susceptibility: 

X = Xtop +XA (5.28) 
Xtop = ^<V^g...|.4t '"^^y^U |^,...> (5.29) 



XA < 



2 r 



NA L 

< %bg,s.\A^\ai >< ai\A\%bg.s. >] (5.30) 



where \ai > stands for the state(s) of the topological multiplet connected 
to the absolute ground-state by A. Using the local properties of A, A^\ai > 
is in the same topological sector as \ai > and < > is at most of 

0(A^x 2-^/2) (see paragraph lS.S.lj) . As Ea^ —Eg.s. is supposed to decrease as 
exp{—L/^) (see Fig. 15. 4|) xtop goes to a constant when the size of the sample 
goes to infinity, provided is small enough ^. In a system with exponentially 
decreasing correlations, decreases a.sl/N and XA is trivially constant at 
the thermodynamic limit. 

In such a phase an infinitesimal field cannot induce a symmetry breaking 
and there could not be any spontaneous symmetry breaking. 



5.6 Other approaches of type I RVB Spin Liquids 

As already said, Sachdev in a large A^, Sp{N) analysis of the Heisenberg 
model on the triangular lattice equally found in the extreme quantum case 
an RVB ground-state and deconfined spinous j95] . 

An interesting mapping of the dimer models on Ising model in transverse 
field has been studied by many authors (See refs. |125| I118j and references 
therein) . This approach is one of the ways to give evidence of the parentage 
between this RVB spin liquid with its fractionalized excitations and the 
deconfining Ising Gauge theories (EEl 1123 II2H1 HIHl EEl 

^Strictly speaking should be < log 2 but Inea. 15.91 is a dramatic overestimate of the 
scalar product in the case of an RVBSL |81| . The reason is that if two dimer coverings 

and c~ maximize (c^|c~) they only differ along a single large (~ L) loop. They have 
different local correlations along the loop and their energy difference is of order L and it 
is very unlikely that their weights in the states i/)"'"(c"'") and {c~) are both of order one. 
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5.7 Summary of the properties of type I RVB Spin 
liquids 

To conclude let us summarize the properties of a type I RVB Spin liquid: 

• It is a phase which does not break either SU(2) symmetry nor any 
spatial symmetries of the lattice. Its ground-state is unique up to 
a topological degeneracy which exists only in systems with an odd 
number of spin-1/2 in the unit cell, living on a 2-torus (more generally 
the degeneracy is 2^, with g the genus of the torus). In that sense it is 
awkward to call such a phase a disordered phase! None of the classical 
ideas associated to disorder are relevant to understand the properties 
of this RVB Spin Liquid phase. If we have to compare it to a liquid 
phase it is more the superfluid phase of '^He that we should have in 
mind! 

• All correlation functions in local observables have only short range 
order, and consequently the susceptibility associated to any local ob- 
servable is zero a T=0. 

• This phase has a gap for all excitations, either in the singlet or the 
triplet sectors and it supports fractionalized excitations (the "spinous" ) . 
The first excitations in the singlet sector correspond in the gauge the- 
ory language to the bosons of the gauge field that Senthil and M. P. 
A. Fisher call "visons" jl27l I128j . Due to the properties of the ex- 
citations we expect them to form continua in the spin sectors. The 
neutron experiment of Coldea and co-workers on CS2CUCI4 is perhaps 
the first experimental proof of such a state |94j . 

• An RVB spin liquid state is expected in presence of competing and 
frustrating interactions. The bandwidth in which this phenomenon can 
be observed is strongly reduced with respect to the original couplings 
and often only a small fraction of the original couplings: the case of the 
MSE system at Jl = —2 and J4 = 1 (the point where all the results 
given here have been calculated) is in some sense rather exceptional, 
as at this point the spin gap is almost as large as the cyclic 4-spin 
exchange. 
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Chapter 6 



Resonating Valence Bond 
Spin Liquid (Type II) 

Type I RVBSL have a gap in the singlet sector. In the Ising gauge theory 
approach this gap is essential for the consistency of the theory |128j . The 
Ising gauge field quasi particles called visons by Senthil and co-workers are 
vortices of the gauge field, they carry a gauge flux, no spin and have long 
range interaction with spinous. If the spectrum of these particles has a gap 
then the spinous are unconfined and the phase is "fractionalized" . If they 
condense, the long range interaction between them and the spinous frustrate 
the motions of the latter which remained confined. The gap in the singlet 
sector (above the topological degeneracy) is a crucial ingredient of Type I 
RVBSL: as we have learned in the previous chapter it is a property of the 
MSE Hamiltonian near the ferromagnetic phase. 

As displayed in Fig. 15. 2[ the MSE Hamiltonian near the antiferromag- 
netic three sublattice Neel order displays a second Spin Liquid phase, with 
no long range order in any observable, but no gap in the sing let sector [Tn7| . 
Such a behavior is not a standard Spin Liquid behavior. Could it be a 
quantum critical behavior similar to that of the Rokhsar Kivelson model at 
y = J? Keeping in mind the restrictions inherent on small size approach 
of such phenomena, we nevertheless think that the observed phenomenon is 
distinct from the critical RK behavior on the square lattice (see below). In 
view of the similarities between the phase seen in the MSE model and the 
phase observed in the Heisenberg model on the kagome lattice, I took the 
step to describe the physics that can be observed at that point as something 
different: a RVB Spin liquid of type II. 

In fact as most of the studies done in many groups have been devoted 
to the Heisenberg model on the kagome lattice (noted HK model in the 
following), it is on this last model that we will center this last chapter. 
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Figure 6.1: The kagome lattice: a lattice of corner sharing triangles. In 
Japanese kagome is a common word which designs both a basket and the 
special canwork for baskets. This is equally a trademark for a popular 
tomato sauce. 

6.1 Miscellaneous models on the kagome lattice 

There has been a large number of studies devoted to different antiferromag- 
netic models on the kagome lattice, which is a 2-dimensional lattice of corner 
sharing triangles (Fig. 16. If) . 

The next neighbor Ising model on such a lattice is disordered, its entropy 
is very large Sj^^^^ = 0.502, more than half the independent spin value, much 

larger than the triangular lattice value S^^-"'^ = 0.323 and of the order of 
Pauling approximation for independent triangles SpauHng = 0.501 |129j . 
This suggests that the correlations in this system are very weak: the model 
remains disordered at all temperatures jl30l I131j . 

Moessner and Sondhi have studied this Ising model in a transverse mag- 
netic field (the simplest way to include in the model some quantum fluctu- 
ations): the model fails to order for any transverse field, at any tempera- 
ture [TH^IT^ . 

The n.n. classical Heisenberg model on the kagome lattice has also an 
extensive T=0 entropy: this is a property easily understood, shared with the 
same model on different lattices with corner sharing units as the checker- 
board lattice or the true three dimensional pyrochlore lattice. On all these 
lattices the n.n. Heisenberg Hamiltonian can be rewritten as the sum of the 
square of the total spin 5*0, of individual units a (a tetrahedron in the 2-d 
and 3-d pyrochlore cases and a triangle for the kagome lattice), which share 
only one vertex. Classical ground-states are obtained whenever Va Sq, = 0. 
This condition fixes the relative position of the three classical spins of a 
triangle at 120 degrees from each other in a plane. But it does not fix 
the relative orientation of the plane of a triad with respect to the planes 
of triads on corner sharing triangles: the model has thus an infinite local 
degeneracy and an extensive entropy jl33lll3l| . This state too is reluctant to 
order: nevertheless thermal fluctuations select coplanar configurations |133[ 
11311 1134j . without long range order in the plane. The order parameter of 
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such a phase is the direction of the local helicity (sometimes called vectorial 
chirality) and defined by its components as: 



where e"^'*' is the antisymmetric tensor. This kind of order is by analogy to 
liquid crystals sometimes called a nematic order. The existence of such an 
order parameter is probably not without relation with the instability of the 
classical Heisenberg model on the kagome lattice to Dzyaloshinsky-Moriya 
interactions ji;-i5| . 

6.2 The next-neighbor spin-1/2 Heisenberg model 
on the kagome lattice: an extreme play-ground 
for "quantum fluctuations" 

The n.n. spin-1/2 quantum Heisenberg model on the kagome lattice has 
equally been the object of many studies 1^1 1?^ TO [TlOl ITITI OHl TO 
ri43. 144 . From these studies one can remember the following facts: 

6.2.1 Ground-state energy per spin 

The Heisenberg model on the kagome lattice has an extremely low energy 
per bond (< Sj.Sj >perbond= —-437) ~ 87% of the energy per bond of 
independent triangles. On this lattice the energy per bond of the spin-1/2 
system is much lower than the classical energy ~ 1.74, a ratio much 

^cl. 

larger than in any other 2-dimensional magnet, that can only be compared 
to the value obtained for the Bethe chain (1.77) (see Table. [LT|) . The kagome 
lattice is the 2-dimensional lattice which offers the larger stabilization due 
to quantum fluctuations. 

6.2.2 Correlations 

Spin-spin correlations jl39j . dimer-dimer correlations (Fig. 16. 2() . chirality- 
chirality correlations jl38j are short range, which is consistent with the pre- 
vious point. 

6.2.3 Spin-gap and absence of gap in the singlet subspace 

There is probably a spin-gap of the order of l/20th of the coupling con- 
stant |2n]- In view of the smallness of this spin-gap with regards to the 
available sizes caution is necessary. The above conclusion was drawn from 
the raw data of exact spectra of samples with up to 36 spins (see Fig. 16.6)1 : 
more precisely from the measurement of Eq{S = Smin + ^) — Eq{S = Smin), 





(6.1) 



on a triangle 
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Figure 6.2: Dimer-dimer correlations in the ground-state of the kagome 
Heisenberg model (blue squares) and in the MSE model (red bullets) versus 
distance. Although the decrease of these correlations is weaker in the KH 
model than in the MSE model, it is nevertheless roughly exponential in these 
first two decades, as the spin-spin correlations are |138| ITr 



where Eq{S) is the lowest energy in the S sector, and Smin = or 1/2 
depending on the parity of the number of spins of the sample. The size 
effects on these results are an order of magnitude smaller than in an ordered 
Neel antiferromagnet. Nevertheless they are still not negligible for these 
sizes. We have thus tried an indirect measurement of the spin-gap along the 
following line. 

The lowest exact eigenstates in each S sector of a sample of N sites define 
the energy per spin of the sample at T = as a function of its magnetization 
m = S/{N/2). For low value of the magnetization one can fit this energy 
per spin to an analytic law of the form ^ : 

e{m) = e{0) + am + bm^/2 + O{m^) (6.2) 

The a and b coefficients depend on N. Their physical significance is clear: a 
measures half the spin gap and 

where x is the homogeneous susceptibility of the medium for fields larger 
than the critical field He {He in convenient units is equal to the spin gap). 

^This analytic form cannot extend beyond m — 1/3, where an angular point appears 
with a discontinuity of the first derivative signaling a magnetization plateau [?. I146| 
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Figure 6.3: Finite size scaling of the spin gap in the spin-1/2 Heisenberg 
model Ti. = Yl<ij> Si-Sj on the kagome lattice. Data are obtained by the 
indirect procedure described in the text. 

This indirect determination of the spin gap leads to a strong renormal- 
ization of the data far small sizes (explaining the large error bars for small 
sizes data in Fig. 16. 3|) . On the other hand the dependence of the renor- 
malized data is now much weaker as can be seen by a comparison of the raw 
results of [23] and Fig. 16.61 with those of Fig. 16.31 A linear extrapolation 
versus 1/N (which should give a lower bound of the spin-gap) leads to the 
value 0.06 for the spin gap (consistent extrapolations give Coo = —0.4365 
and Xoo = 0.34). All these determinations are in agreement with the re- 
sults obtained from the raw data. Nevertheless is should be underlined that 
even for N=36 at the smallest non zero magnetization the linear term of 
equation (|6.2)) . is only 90% of the quadratic term: this shows the limit of 
confidence in our assumption on the existence of a spin-gap ^. 

6.2.4 An exceptional density of low lying excitations in the 
singlet sector 

Whatever the ultimate fate of the spin gap a still larger surprise emerges 
from the exact spectra: the absence of gap in the singlet sector and the 
anomalous density of low energy states adjacent to the ground-state. We 
have measured the number of singlet levels in the spin-gap (taken as a natu- 
ral energy band- width of the problem): this number increases exponentially 
fast with N as 1.15^ (see Fig.lfTl). 

The first immediate consequence of this property is the existence of a 
T = residual entropy in this model. This came as a shock for many 
scientists who had the idea that the quantum dynamics "should" lift the 
degeneracy of the incipient Ising model or dimer models (see next section). 

^This spin gap result is in my point of view the less reliable result among all those 
described in these lectures. In view of the spin-spin correlations which seems exponentially 
decreasing and not critical |189) . and of the fact that results between 27 and 36 may signal 
a cross-over behavior in the finite size effects, we will not question this point further for 
the moment. It would nevertheless be useful to develop a very precise analysis of finite 
size effects between 27 and 36 (as we have done for the MSE model see Fig. I5.3|l to try to 
confirm this conclusion. 
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N 

Figure 6.4: Logarithm of the number of singlet states in the spin gap ver- 
sus sample size (black squares). The short dashed and long dashed curves 
display the theoretical law (Eq. 16. 5|) (short dashes: p = 1, long dashes: 
p = 2). 

This does not seem to be the case!! 

Some remarks are necessary to fully appreciate this property. 

The total number of states of a sample of N spins 1/2 is 2^. These states 
are stretched on an energy scale of the order of NJ where J is the coupling 
constant of the Hamiltonian. This implies that on most of the spectrum 
the density of states increases exponentially with N. If we specialize to the 
5 = sector as we will do below, the picture is not very different: the 
number of states is — ~ ^(iv") i^e^e too, in most of the 

spectrum the density is exponentially increasing with N. 

BUT in all the phases that we have studied up to now, the nature of the 
ground-state and of the low lying excitations leads to a different behavior 
at the edges of the band. Typically the ground-state degeneracy is 0(1) in 
VBC and in type I RVBSL and it is 0{N°') in Neel ordered states with a 
sublattices. In all these situations the low lying excitations can be described 
as modes or quasi-particles. In the corresponding energy range one typically 
counts levels associated with one quasi particle excitations. This always 
leads to density of states increasing as a power law as a function of N. 
Inclusion of multi-particle excitations can be done in an average way: let us 
suppose that single particle excitations have a dispersion law: 

uj{k) = F. (6.4) 

In a d-dimensional space the internal energy of such a system increases with 
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temperature as T^^^'^)/^, the specific heat as T'^/p as the entropy S. For a 
system with A'^ spins the average energy range W excited at the temperature 
T is of the form: W = C'^^NT^^^'^^^'p . This imphes that the entropy depends 

d 

on A'^ as: S (x N (■^) ''+'' and thus the average number of excited levels 
increases as: 

ln{Af) oc (6.5) 

3 

As an example, let us consider the quantum critical phase of the Rokhsar 
Kivelson model on the square lattice at the quantum critical point (Chapter 
4 page 60): in that case the critical correlations decrease as r~^, the disper- 
sion law is linear in k, the logarithm of the number of states increases as 
jyi/s (short dashes of Fig. 16. 4j) . much slower than the exact results. 

Even in including many particle excitations one would thus expect a 
number of levels increasing much more slowly than in the exact results! 
Except if we accept that there is infinitely soft low energy modes {p oo), 
then we recover the correct density of low lying levels and a residual entropy. 
It is still unclear if we can do a connection between the "zero mode" of the 
classical model at T=0 and this picture. And we cannot completely indulge 
ourselves in saying that quantum fluctuations are unable to lift the classical 
degeneracy as quantum fluctuations seems to open a spin gap! 

A physical consequence of this exceptional density of low lying singlets 
can be observed in the specific heat (Fig. 16. 5() : at low temperature the 
specific heat of this spin system is unusually insensitive to large magnetic 
fields. This is easily understood if we suppose that in this energy range 
the excitations are essentially singlets |143j . This result is to be compared 
to the experimental results of Ramirez et al. jl4Sj on SrCrGaO where the 
specific heat around 5K has an extremely low sensitivity to magnetic fields 
up to 10 Tesla, whereas the homogeneous susceptibility in this range of 
temperature is probably very low if we notice that it turns down around 

50K mni. 

6.2.5 Anomalous density of states in other spin sectors 

This anomalous density of low lying states has equally been observed in 
the spin 1/2 sector (where the law could be fitted to A^1.15^), in the spin 
1 sector as well as in other sectors with larger total spin. It should be 
noticed that such a density of states implies the absence of an intrinsic 
energy scale for the low lying excitations: a phenomenon that has been 
observed in inelastic neutron scattering (ref. |15nj and refs. therein) and 
theoretically in the imaginary part of the dynamic susceptibility calculated 

thank S. Kivelson for giving me this idea for the computation of the muhi-particle 
density of states. The errors, if any, are mine. 
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within the dynamical mean field theory |151j . A very high spin susceptibility 
just above the spin gap is not excluded in spin-1/2 compounds jl52j . 



6.2.6 Unconfined Spinons 



The spinon binding energy (Fig. I6.6j) has been computed in the HK model 
along the same line as it has been done in Chapter 5 page 75 for the MSE 
model. The only difference lays in smaller and smoother shape effects in the 
HK case, which allows to do analytic fit of the Gncrgy per spin 6(*S''mm5-^) 
versus even and odd samples separately. Eq. 15. 6[ is then written, 

with reference to the interpolated value of the energy of even samples for odd 
value of N: this takes into account the finite size effects with more subtlety 
than the (2A^ + l)eoo term of equation 15.61 All these features explain the 
smoother behavior of the spinon binding energy shown in Fig. 16.61 From this 
data one can rather safely conclude that spinons are probably unconfined in 
the HK model. 

The global picture of this phase is thus that of a Spin Liquid, with 
no long range correlations in local observable, unconfined spinons and a 
residual entropy of singlets at T = 0, which is one of the manifestations of 
an extraordinary large density of states in each S subspace. 
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Figure 6.5: Specific heat (dot-dashed curve), entropy (dotted curve) and 
spin-susceptibility (full line) of the Heisenberg model on the kagome lattice 
(exact diagonalization on a N=18 sample, and Padde approximants to high 
temperature series jl47| ). The third bump in the specific heat at very low 
temperature is an artifact of the small size and the associated discretization 
in the singlet sector. We think that the intermediate bump is a real feature 
of the specific heat which subsists up to the thermodynamic limit |143j . 
Note the entropy of the singlets in a range of temperature below the spin 
gap where the spin excitations are negligible as the spin-susceptibility. 
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Heisenberg model on the kagome lattice: 
Spin Gap and Spinon binding energy 
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Figure 6.6: Finite size scaling on the spin gap: red squares (raw data) and 
the spinon binding energy: blue squares (unpublished results). 

6.3 Next-neighbor Resonating Valence Bond de- 
scription of the spin-1/2 kagome antiferromag- 
net 

Considering a supposed-to-be large spin-gap, Zeng and Elser jl4Uj proposed 
a description of the ground-state and low lying excitations of the kagome 
model in the basis of next neighbor Valence Bonds. We know now that 
the spin gap is certainly smaller than it was expected in 1995, nevertheless 
Mila and Mambrini jl42l I144j have convincingly shown that the picture 
of a next neighbor resonating Valence Bond Spin liquid captures some of 
the most perplexing features of this magnet and specifically the absence 
of gap in the singlet sector and the exponential number of singlets in any 
given range of energy. This probably implies the absence of an intrinsic low 
energy scale, which is consistent with the thermal behavior of the dynamic 
spin susceptibility calculated by Georges and coll. jl51j . This feature is 
typical of a critical state, but as seen in the above discussion, the simple RK 
picture does not seem to fit nicely to the exact diagonalization data: may 
be the available sizes are too small or the behavior of this quantum system 
corresponds to something definitely new. 

More generally this picture of the ground-state and first excitations as 
resonances between an exponential number of dimer coverings gives a qual- 
itatively interesting picture of the low temperature physics of different ox- 
ides that can be described as kagome antiferromagnets. The low temper- 
ature specific heat of SrCr9Gai20i9 is apparently dominated by local sin- 
glet states |148) . The magnetic excitations of this same compound as seen 
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by muons can be described as spins 1/2 itinerant in a sea of singlets |25j . 
SrCr9Gai20i9 exhibits at about 5 K a very large increase of its non linear 
susceptibility reminiscent of spin glasses (153] , but neutrons and muons show 
that most of the spins are not frozen below this temperature and exhibit still 
very rapid fluctuations |154j . The same phenomena have been observed in 
two jarosites that are equally good models of kagome antiferromagnets with 
half-odd- integer spin per unit cell |155l I156j . 

6.4 Haldane's conjecture 

Whereas the classical Heisenberg model on the kagome, checker-board and 
pyrochlore lattices share the property of an extensive entropy and disorder 
at T = 0, their quantum counterparts are quite different. As it has been 
explained in Chapter 4, the Heisenberg model on the checker-board lattice 
has an ordered Valence Bond Crystal with gaps to all excitations. 

Less is known on the ground-state of the Heisenberg model on the 3- 
dimensional pyrochlore lattice: Canals and Lacroix jl57j have shown that 
their spin-spin correlations are short ranged and they have seen on a 16 
sites spectrum that the first excitations were singlet ones. Having done the 
spectrum of a 32 sites pyrochlore sample, we confirm that the first excitations 
are still singlets for this size. But the finite size effects on these excitations 
between 16 and 32 are very large and it remains possible that the gap in the 
singlet sector be larger than the singlet-triplet gap in the thermodynamic 
limit [?, ?]: in any case there is definitely a gap in the singlet sector in this 
last model!'^ 

On the other hand, recent results from Hida jl6()j show that there is a 
gap to all excitations in the S=l HK model. 

All these results seem to confirm Haldane's conjecture: among these 
strongly frustrated systems with an extensive degeneracy in the classical 
limit the spin- 1/2 kagome antiferromagnet is the only system to have an 
half- odd- integer spin in the unit cell and gapless excitations. The spin-1/2 
Heisenberg model on the checkerboard lattice or on the pyrochlore lattice 
and the spin-1 Heisenberg model on the kagome lattice have integer spins 
in the unit cell and quantum fluctuations lead to gapful excitations. 



^Tsutenegu has recently developed an effective description of the singlet sector, where 
he has a soft mode in the singlet sector. I am a bit skeptical on his approach which 
strongly and artificially breaks the symmetries at the mean-field level |158| and then treats 
in a semi-classical approach the fiuctuations |159) . The Core approach of E. Berg and 
collaborators [?] seems more appropriate to deal with these systems where the quantum 
dimerization is probably the dominant phenomenon 
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